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ABSTRACT 


Difficulties  arising  from  inherent  inaccuracies  in  structural 
modelling  and  from  the  high  dimensionality  of  the  dynamical  system 
force  a  re-examination  of  the  problem  of  optimal  control  of  large 
flexible  structures  within  the  context  of  stochastic  system  theory 
In  this  report/  the  design  of  active  structural  control  is  formu¬ 
lated  as  the  mean-square  optimal  control  of  a  linear  mechanical 
system  with  stochastic  parameters.  In  practice,  a  complete  prob¬ 
abilistic  description  of  modal  parameters  can  never  be  provided, 
and  a  suitable  design  approach  must  accept  very  limited  a  priori 
data  on  parameter  statistics.  In  consequence,  we  formulate  the 
mean  square  optimization  problem  using  a  complete  probability 
assignment  induced  by  the  available  data  through  use  of  a  maximum 
entropy  principle.  Furthermore,  it  is  proposed  to  acknowledge  as 
available  the  minimum  set  of  a  priori  statistical  data  on  parame¬ 
ter  variations  which  is  needed  to  preserve  any  measure  of  model¬ 
ling  fidelity. 

To  fix  ideas,  we  specifically  address  the  problem  of  full 
state  feedback  regulation  of  a  linear  structural  system  with 
statistical  variation  only  in  the  open  loop  frequencies.  Examin¬ 
ing  the  phenomenology  of  modal  frequency  uncertainties  we  discern 
the  so-called  "modal  decorrelation  times"  as  the  minimum  data 
required  to  preserve  2nd  moment  response  characteristics  at  high 
levels  of  uncertainty  or  for  high  order  structural  modes. 

The  decorrelation  times  are  closely  related  to  the  damping 
times  of  the  parameter  ensemble  averaged  modal  response  and  their 
reciprocals  constitute  fundamental,  albeit  unconventional, 
measures  of  the  variation  of  the  modal  frequencies  about  their 
nominal  or  mean  values.  Choosing  to  acknowledge  the  decorrelation 
times  as  the  available  data,  the  complete  probability  assignment 
which  is  otherwise  maximally  unconstrained  is  a  white  parameter 
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model  in  which  the  noise  intensities  are  inversely  proportional 
to  the  decorrelation  times.  The  mean-square  optimization  problem 
therefore  reduces  to  the  solution  of  a  stochastic  Riccati  equation 
of  a  form  arising  from  the  state  dependent  noise  problem. 

Certain  features  of  the  stochastic  Riccati  equation  are  next 
explored.  It  is  shown  that  under  weak  restrictions  a  unique 
positive  semi-definite  solution  exists  for  all  values  of  the 
decorrelation  times.  Also,  2nd  moment  stability  is  guaranteed 
for  the  closed- loop  system.  Thus,  the  need  for  design  iteration 
to  ensure  robustness  with  respect  to  stability  is  largely  elimi¬ 
nated  within  the  proposed  approach. 

Furthermore,  in  the  limit  as  all  decorrelation  times  approach 
zero  (i.e.,  uncertainties  in  all  modal  frequencies  increase  with¬ 
out  bound) ,  the  solution  of  the  stochastic  Riccati  equation  yields 
a  rate  feedback  control  law  which  is  stable  for  all  values  of 
modal  frequencies  or  damping  ratios.  In  the  more  typical  case  in 
which  low  order  structural  modes  are  relatively  well  known  while 
modelling  accuracy  degrades  for  higher  order  modes,  the  stochastic 
Riccati  equation  produces  a  control  which  approaches  the  asymp¬ 
totic  rate  feedback  form  for  high  order  modes  and  yet  closely 
resembles  "high  authority"  deterministic  plant  design  for  the  low 
order,  relatively  well-known  modes.  These  two  regimes  exist  as 
limiting  qualitative  features  of  a  global  control  law  for  which 
stochastic  stability  is  guaranteed. 

The  above  qualitative  features  permit  a  numerical  scheme  for 
determination  of  the  optimal  gain  matrix  in  which  the  computational 
burden  is  mainly  associated  with  the  relatively  well-known  or 
"coherent"  modes.  As  long  as  the  "coherent  system"  is  of  modest 
dimension,  the  stochastic  Riccati  equation  admits  of  practical 
numerical  treatment  for  systems  of  arbitrarily  large  order.  Thus, 
the  proposed  design  approach  eliminates  the  need  for  modal  trunca¬ 
tion  with  its  attendant  spillover  instability  problems. 
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OPTIMAL  REGULATION  OF  STRUCTURAL  SYSTEMS 
WITH  UNCERTAIN  PARAMETERS 


1 .  INTRODUCTION 

The  advent  of  the  Space  Shuttle  has  prompted  considerable 
attention to  the  design  and  control  of  large,  lightweight 
space  platforms.  Feasibility  studies  have  identified  very 
stringent  tolerances  on  figure  and  precision  pointing  require¬ 
ments.  (2,3,4) 

In  the  face  of  severe  mission  requirements,  vibration  sup¬ 
pression  by  purely  passive  structural  damping  may  be  inadequate 
and  active  electronic  control  of  the  structure  must  be  contem¬ 
plated.  Due  to  its  significant  influence  on  mission  performance, 
vehicle  elasticity  must  be  carefully  accounted  for,  so  that  many 
degrees  of  freedom  are  required  in  the  system  model.  The  high 
dimensionality  of  the  dynamic  system  renders  classical  control 
design  techniques  excessively  cumbersome.  Thus,  the  more  sys¬ 
tematic  approach  of  "modern"  optimal  control  theory  is  preferable. 

Application  of  modern  control  theory  to  the  design  of  active 
structural  control  faces  two  fundamental  problems.  In  the  first 
place,  control  effectiveness  provided  by  the  optimum  control¬ 
ler  and  estimator  gains  is  sensitive  to  errors  in  modelling  plant 
dynamics,  i.e. ,  slight  variations  in  the  model  parameters  and 
errors  in  environment  specification.  This  is  the  problem  of 
"robustness."  Secondly,  there  is  the  difficulty  posed  by  the 
large  dimension  of  the  dyneimical  model  that  must  be  employed. 

In  particular,  stability  of  the  closed-loop  system  cannot  be 
guaranteed  if  dynamically  significant  modes  are  truncated  from 
the  plant  model.  This  is  the  well  known  problem  of  control  and 
observation  "spillover" . 
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If  the  computational  difficulties  were  not  prohibitive,  one 
might  attempt  to  circumvent  the  above  problems  by  use  of  a  very 
accurate  high-order  plant  model  in  a  standard  Linear-Quadratic 
formulation.  The  analysis  of  high-order  structural  modes,  how- 
ever,  is  plagued  by  two  essential  difficulties.'  '  '  First, 
truncation  of  the  infinite  dimensional  system  to  a  finite¬ 
dimensional  system  is  always  inherent  to  the  finite-element 
representation.  The  corresponding  errors  incurred  tend  to 
accumulate  with  increasing  frequency  and  the  calculated  mode 
shapes  and  frequencies  are  known  to  be  inaccurate  for  higher- 
order  modes,  even  for  idealized  systems.  Secondly,  and  in 
addition  to  mathematical  modelling  difficulties,  the  higher  modal 
frequencies  and  mode  shapes  are  very  sensitive  to  small  details 
of  geometry,  construction  and  material  properties.  Thus,  in 
considering  structural  response  at  sufficiently  high  frequency, 
the  relevant  structural  details  may  never  be  determined  and 
modelled  with  acceptable  accuracy. 

Note  that  if  the  system  order  did  not  pose  a  problem  one 
could  deal  with  parameter  errors  within  a  traditional  design 
philosophy.  First  design  the  controller  by  some  method,  then 
characterize  robustness  with  respect  to  various  system  properties 
by  determining  acceptable  bounds  on  parameter  variations.  Recent 
results '  '  on  the  multivariable  robustness  issue  greatly  facili¬ 
tate  this  process.  Unfortunately,  since  the  number  of  system 
parameters  whose  variation  must  be  considered  increases  rapidly 
with  the  system  order  this  approach  would  seem  to  entail  great 
complexity. 

The  difficulties  posed  by  dimensionality  and  parameter  sen¬ 
sitivity  stem,  in  part,  from  a  reliance  upon  design  methods  which 
implicitly  assume  complete  information  on  system  parameters.  It 
is  felt  that  approaches  which  include  measures  of  parameter  uncer¬ 
tainty  as  part  of  the  control  formulation  offer  significant  theo¬ 
retical  advantage. 


A  number  of  such  controller  design  methods  have  been  advanced 
in  the  last  several  years,  and  comparative  assessments  are  to  be 
found  in  References  (10)  and  (11) .  A  distinct  philosophy  with 
regard  to  the  modelling  of  parameter  uncertainties  is  exemplified 
by  the  guaranteed  cost  control  method  '  which  attempts  to 
bound  the  effects  of  uncertainties.  Here  the  system  parameters 
are  assumed  to  lie  in  a  closed  bounded  region  and  a  modified 
Riccati  equation  is  devised  such  that  closed-loop  system  behavior 
is  acceptable  for  all  value j  of  the  parameters  within  specified 
limits.  Recently  this  has  been  noted  to  produce  large  controller 
gains  and  relatively  large  control  effort  with  overdamped  domi¬ 
nant  closed-loop  poles.  Although  a  special  iterative  procedure 

(13 ) 

has  been  devised  to  remedy  this  disadvantage,  the  applicabil¬ 

ity  of  this  approach  to  systems  having  a  large  number  of  uncertain 
parameters  has  not  been  demonstrated. 

In  any  case,  methods  which  assume  parameter  variations  within 
a  limited  range  are  basically  non-statistical  in  character.  A 
more  thorough-going  approach  would  employ  the  general  concepts  of 
stochastic  optimal  control , allowing  parameter  probability 
distributions  of  unbounded  support.  We  refer  specifically  to  the 
concepts  of  dual  control  and  related  ideas  as  clarified  by 
Bar-Shalom  and  Tse.  Within  this  general  context,  parameter 

uncertainties  are  not  necessarily  accepted  at  their  a  priori 
values;  instead,  a  control  is  considered  truely  optimal  when  it 
has  a  dual  effect  -  i.e.,  when  it  not  only  affects  the  state  of 
the  system  but  also  reduces  the  uncertainty  of  the  state. 

Unfortunately,  the  formulation  of  adaptive  control  is  a  non- 

f  1 8 ) 

linear  stochastic  control  problem'  '  and  approximations  are 
usually  needed  to  achieve  a  practical  solution.  One  attractive 
approximation  replaces  the  actual  system  parameters  by  parameters 
which  are  uncorrelated  in  time,  thereby  making  identification 
impossible  and  precluding  probing  action.  Not  only  is  the  problem 
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rendered  solvable  but  this  approach  is  design  conservative  in  the 

sense  that  it  bounds  the  deterioration  in  performance  due  to 

(19) 

unknown  parameters  and  allows  us  to  obtain  an  inherent  caution 

in  the  control. 

Within  the  non-dual  approximation  and  in  addition  to  the  well- 

(21  22 ) 

known  state  dependent  noise  formulation  for  continuous  systems,  ' 

much  recent  work  has  been  directed  to  the  discrete-time  case.  In 

particular,  the  uncertainty  threshold  principle  of  Athans  and 
(19  23) 

Ku  '  and  recent  results  for  the  general  multivaria’'le  prob- 
( 24 ) 

lem  '  provide  valuable  qualitative  insj.nht  in  the  present  prob¬ 
lem.  Clearly,  the  possibility  of  an  uncertainty  threshold  is  a 

critically  important  issue  for  active  vibration  suppression. 

( 24 ) 

Moreover,  experience  with  the  multivariable  case  shows  that 

the  presence  of  a  large  number  of  uncertain  parameters  can  be  the 
source  of  prohibitive  complexity. 

The  view  espoused  here  is  that  a  suitable  design  method  must 
not  only  be  consistent  with  a  non-dual  approximation  but  must 
incorporate  very  limited  information  on  parameter  statistics 
thereby  reducing  the  number  of  measures  of  parameter  uncertainty 
to  a  manageable  level. 

Actually,  in  the  specific  area  of  response  estimation,  the 
issue  raised  by  dimensionality  and  parameter  uncertainty  and  the 
need  for  a  tractable  description  of  parameter  statistics  have 
been  addressed,  in  part,  by  Statistical  Energy  Analysis.  Moti¬ 
vated  by  a  concern  with  high  frequency  vibration  and  acoustical- 
structural  interaction,  the  development  of  SEA  has  been  largely 
concurrent  with  that  of  modern  optimal  control  theory  and  has 
given  rise  to  an  extensive  literature  (see  References  25  and  26 
for  a  general  review  of  SEA  procedures).  Broadly  speaking,  SEA 
attempts  to  formulate  a  "contracted"  description  of  system 
response  consistent  with  the  inherent  uncertainties  of 
structural-acoustical  modelling. 
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In  its  simplest  form,  SEA  divides  a  complex  structure  into 
"subsystems"  which  are  considered  as  repositories  of  vibrational 
energy.  Each  subsystem  consists  of  a  group  of  "similar"  energy 
storage  modes.  Modes  which  play  a  significant  role  in  transmis¬ 
sion,  dissipation  and  storage  of  energy  and  which  have  nearly 
equal  excitation,  damping  and  coupling  to  other  subsystems  are 
assigned  to  a  particular  group.  Generally,  subsystems  are 
associated  with  a  particular  frequency  band  and  may  be  associated 
with  separate  structural  components.  SEA  achieves  substantial 
simplification  by  using  the  approximate  result  that  power 
flow  between  pairs  of  coupled  modes  is  proportional  to  the  dif- 
ference  in  the  average  modal  energies,^  '  '  the  constant  of 

proportionality  being  termed  the  "coupling  loss  factor."  At  this 
juncture,  the  effect  of  randomness  in  the  modal  frequencies  is 
accounted  for  by  averaging  the  coupling  loss  factors  over  appro¬ 
priately  defined  statistics.  In  taking  this  most  important  step, 
one  typically  renounces  all  knowledge  of  the  relative  modal 
frequency  locations  within  each  subsystem,  and  assumes  these  to 
be  randomly  distributed  over  the  subsystem  frequency  band  (the 
assumption  of  modal  disorder).  In  this  manner,  the  SEA  model 
incorporates  a  degree  of  information  regarding  system  parameters 
which  is  commensurate  with  the  limited  information  actually 
available.  Under  these  assumptions  the  over-all  subsystem  power 
balance  relations  are  the  governing  equations  and  these  involve 
only  the  average  modal  energies  of  the  modes  within  each  sub¬ 
system.  The  average  modal  energies  thus  constitute  the  basic 
dynamical  variables  and  all  other  measures  of  second-moment 
response  may  be  deduced  from  them.  At  high  frequency,  where 
subsystems  contain  many  modes  and  the  uncertainties  in  modal 
frequencies  are  large,  the  use  of  energy  and  power-flow  variables 
results  in  a  drastic  decrease  in  the  number  of  measures  required 
to  characterize  the  response  with  an  accuracy  consistent  with  the 
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uncertainties  in  modal  parcimeters.  Most  importantly,  if  suffi¬ 
cient  modal  density  exists,  the  total  average  energy  of  all  modes 
above  a  given  frequency  may  be  estimated  via  simple  asymptotic 
formulae  and  the  problem  of  dimensionality  is  circumvented. 

Although  SEA  displays  many  advantages  in  estimating  high 
frequency  response,  its  use  may  be  inappropriate  at  low  frequen¬ 
cies.  First,  the  confidence  intervals  associated  with  the  cal¬ 
culated  modal  energies  are  inversely  proportional  to  the  number 
of  modes  within  the  subsystems  or,  equivalently,  directly  pro¬ 
portional  to  modal  density.  At  low  frequency,  the  system  lacks 
sufficient  modal  density  to  allow  predictions  with  an 
acceptable  degree  of  certainty.  Secondly,  the  assumption  of 
modal  disorder,  although  it  produces  a  simple  model,  is  unsuit¬ 
able  at  low  frequency.  This  assumption  would  effectively  discard 
all  the  detailed  and  relatively  accurate  information  describing 
the  low  frequency  modes.  On  the  other  hand,  such  information  is 
of  great  advantage  in  controller  design  and  must  not  be  ignored. 

Although  the  specific  procedures  of  SEA  may  not  be  directly 
applicable,  the  underlying  philosophy  offers  an  attractive 
approach  to  the  present  problem.  In  essence,  the  SEA  approach 
to  response  estimation  is  to  incorporate  incomplete  system  infor¬ 
mation  within  the  dynamical  model  by  limiting  consideration  to 
randomness  in  modal  frequencies.  Since  it  deals  with  performance 
measures  defined  over  the  entire  system  ensemble,  SEA  incorporates 
the  effect  of  parameter  uncertainties  at  a  fundamental  level.  But 
most  importantly,  this  modeling  approach  has  profound  consequences 
for  the  problem  of  dimensionality.  In  optimal  regulation  as  in 
response  estimation,  the  "curse  of  dimensionality"  is  manifested 
in  the  great  mass  of  processing  of  fundamental  data  (of  system 
models  presuming  complete  information)  required  for  response  cal¬ 
culation  or  formulation  of  an  optimal  control  policy.  The  example 
of  SEA  intimates  the  possibility  that  by  use  of  models  which 
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include  limited  system  information,  we  may  so  arrange  matters 
that  the  processing  required  for  control  policy  formulation  may 
be  similarly  limited. 

Drawing  inspiration  from  the  essential  ideas  of  SEA,  this 
report  sets  forth  the  basis  of  a  formulation  which,  it  is  hoped, 
will  circumvent  the  difficulties  faced  by  current  design  methods. 
Section  2.1  formulates  the  problem  of  continuous  time,  full 
state-feedback,  linear  optimal  regulation  of  a  structural  system 
with  uncertain  modal  frequencies.  Next,  it  is  recognized  that 
the  statistical  model  of  parameter  uncertainties  must  be  derived 
from  severely  limited  data.  In  Section  2.2  we  identify  a  minimal 
data  set  (the  modal  "decorrelation  times")  which  preserves  certain 
asymptotic  properties  of  the  open-loop  system  in  the  case  of  large 
uncertainties  or  for  high  order  modes.  Choosing  to  acknowledge 
only  this  data  set,  the  statistical  model  which  is  otherwise  max¬ 
imally  unconstrained  is  determined  by  use  of  an  entropy  principle. 
In  consequence  we  obtain  a  white  parameter  model  depending  on 
relatively  few  measures  of  frequency  uncertainty.  Section  2.3 
sets  forth  the  solution  for  optimal  linear  regulation  and  a 
special  form  of  the  stochastic  Riccati  equation  arising  from  the 
state-dependent  noise  problem  is  obtained.  The  section  concludes 
by  summarizing  the  overall  rationale  of  the  proposed  design 
approach. 

Section  3  examines  the  properties  of  solutions  to  the 
stochastic  Riccati  equation.  Specifically  addressed  are  the 
issues  of  closed-loop  stochastic  stability  and  the  influence  of 
uncertainties  on  the  effective  dimensionality  of  the  system.  In 
particular,  it  is  found  that  the  computational  burden  may  be 
reduced  to  that  associated  with  the  relatively  few  "well-known" 
modes,  permitting  the  use  of  high  order  models  and  largely  elimi¬ 
nating  the  need  for  modal  truncation. 


2.  STOCHASTIC  OPTIMAL  CONTROL  UNDER  LIMITED  PARAMETER 

INFORMATION 

2 . 1  Problem  Statement 

To  begin  formulation  of  the  linear  regulator  problem  it 
is  most  suitable  to  write  the  equations  of  motion  for  the  struc¬ 
ture  normal  mode  coordinates.  The  state-space  form  of  these 
equations  may  be  written: 


X  =  AX  +  Bu  +  w  ,  t  e[tQ»tj^] 


x(t^)  =  0 


A  e  B  e  u  e  R^ 


(1) 


where,  for  convenience,  we  restrict  consideration  to  zero  initial 
state.  Here  x  is  the  vector  of  modal  coordinates  and  velocities 
with  its  odd  indexed  elements  representing  modal  displacements 
and  the  adjacent  even  indexed  elements  giving  the  corresponding 
modal  velocities.  With  this  convention  and  supposing  that  linear 
gyroscopic  terms  may  be  neglected  and  that  damping  is  propor¬ 
tional,  the  dynamic  matrix  A  assumes  the  form: 


A  =  block-diag 
k=l , . . . ,n 


0 


1 


-2nj,a)k 


(2) 


where  and  rij^  are  the  k^^  modal  frequency  and  damping  ratio. 
To  simplify  the  developments  of  this  section  and  without  loss 
of  generality  we  consider  only  the  elastic  modes  of  the  system 
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so  that  all  the  are  non-zero.  Henceforth,  we  assume  small  but 
non-zero  structural  damping,  i.e. , 

0  <  <<  1 

for  all  k. 

Furthermore,  all  the  odd  numbered  rows  of  the  input  matrix  B 
are  zero: 


where  the  non- zero  elements  are  proportional  to  the  normal  mode 
shapes  at  2,  actuator  locations. 

w  is  a  vector  of  white  noise  disturbance  forces.  Because 
of  the  convention  regarding  our  state  space  representation,  the 
intensity  matrix,  v,  assumes  the  form: 

\j  =  0  '■  3  odd  (4) 

where  this  is  assumed  non-negative  definite. 

The  control  input  vector,  u,  is  assumed  to  arise  from  linear, 
full  state  feedback; 

u  =  -K(t)  X  (t)  (5) 

where  K(t)  is  a  time  varying  gain  matrix. 


A.  . 


( 29 ) 

The  standard  linear  quadratic  formulation,  assumes  com¬ 

plete  and  accurate  knowledge  of  the  system  dynamics  (A,  B,  V 
precisely  known)  and  seeks  the  K(t)  which  minimizes  the  quadratic 
functional 


J  k  E[/^  dt[x'^Rj^x  +  u'^R2u]  ] 
^o 

Rj^  >  0  ,  R2  >  0 


a. 


>> 


b.  J 


(6) 


subject  to  the  constraints  imposed  by  (1) .  To  fix  ideas  we 
have  not  included  terminal  state  weighting  since  this  involves 
only  trivial  modification  of  the  following  results. 

As  a  further  preliminary  step,  it  is  convenient  to  express 
the  above  relations  in  the  eigen-basis  of  the  uncontrolled  (K=0) 
system.  In  view  of  the  assumption  of  small  damping  we  may  sim¬ 
plify  this  process  by  introducing  the  resonant  approximation 
for  A: 

A  =  block-diag 
k— 1,. . . ,n 


•^k'^k 


-u). 


■^k‘"k 


(7) 


so-called  because  the  difference  between  damped  and  undamped 
natural  frequencies  is  neglected.  The  eigenvalues  of  (7)  differ 
from  those  of  (2)  by  terms  of  second  order  in  the  damping,  while 
the  eigenvectors  differ  by  terms  of  first  order.  Within  this 
replacement,  the  eigenvector  matrix  of  A  is 

4>  =  block-diag 
k=l , . . . ,n 


ia,k 


(8) 


L 
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A. 


Def ining : 

=  <t>  X  ,  W  =z  <P  w 


(9, 


and 


u  ^ 

$  ^A<li  =  diagiaij^  (i-Hj^)  ,Wj{-i-nj^) 

a. 

...u,n(i-n„),u.n(-i-n^)} 

W 

11  L> 

b. 

► 

^  A 

K4> 

c. 

< 

II  D> 

.-1„.-1H 

9  V'f 

d. 

.,5 

e-  - 

the  optimal  linear  regulator  problem,  equations  (1)  through  (6) , 
may  be  re-stated  as: 


min : 

< 

J 

=  E(/ ^  dt  f“[o.  + 

^O  L  ^ 

a. 

K  real 

0^  ^  0  ,  >  0 

b. 

• 

f; 

=  (u-e<)c  +  w  ,  t 

c . 

o^ 

II 

o 

d. 

where  w  is  white  noise  with  intensity  matrix  v. 

^  Now,  the  inevitable  errors  in  the  finite  element  model  induce 
statistical  variation  in  all  the  matrices  defined  in  (10) .  How¬ 
ever,  to  gain  an  analytical  foothold  on  the  problem,  we  begin  by 
consideration  of  the  simplest  possible  case  in  which  only  the 
open-loop  eigenvalues  (the  u  matrix)  are  subject  to  random  varia¬ 
tion.  Note  that  whatever  the  modelling  errors,  the  real  parts  of 
u  will  be  negative  and  small.  Thus  it  is  likely  that  the  impact 
of  random  fluctuation  in  the  open-loop  frequencies  on  closed-loop 
performance  will  overshadow  randomness  in  the  damping  terms. 
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Consequently,  we  further  limit  consideration  to  random  Im(y). 
Specifically,  it  is  assumed: 


u  =  y  +  v(t) 


a. 


U  =  diag{a)j^  (i-Hj^)  (-i-Hj^)  , —  ^ 

v(t)  4  diag  {i  Im(u.)d.  (t)} 

k=l,...,2n 


b. 

c. 


>  (12) 


where  the  uij^,  k=l.  .  .n  are  the  nominal  or  mean  values  of  the  modal 
frequencies.  The  6j^(t);  k=1...2n  are  assumed  real  valued,  zero- 
mean  and  stationary  random  processes  in  time  and  are  mutually 
statistically  independent  and  independent  of  the  disturbance 
noise  w. 

The  above  model  is  consistent  with  an  SEA  point  of  view. 

The  standard  SEA  formulation  tacitly  acknowledges  that  the 
energetics  of  modal  interactions  are  most  sensitive  to  relative 
frequency  locations  so  that  randomness  in  the  frequencies  is 
usually  emphasized. 

Obviously,  a  general  treatment  would  require  that  all  matri¬ 
ces  be  random.  However,  the  above  restrictions  offer  an  appro¬ 
priately  simple  point  of  departure  and  permit  relatively  easy 
interpretation.  Moreover,  as  will  ultimately  be  seen,  the  very 
special  problem  considered  here  still  exhibits  important  features. 

To  include  the  effect  of  frequency  uncertainties  while  retain¬ 
ing  the  form  of  linear  quadratic  theory,  we  again  employ  (11. a) 
but  extend  the  averaging  operation  over  the  parameter  ensemble  as 
well  as  over  the  disturbance  ensemble.  Having  done  this,  it  is 
advantageous  to  re-phrase  (11)  in  terms  of  the  co-state  matrix, 

0,  of  the  associated  deterministic  plant  regulation  problem. 

Then,  within  the  above  restrictions,  the  problem  is  to  determine 
a  K  (such  that  K  is  real)  to  minimize  : 


subject  to  the  constraints 

P  =  Efp] 

-p  =  (u  +  v(t)-g<)^p  +  p  (p  +  v(t)-e<) 

+  +  <  R2< 

P(t^)  =  0 

That  this  statement  follows  from  (10)  and  (11)  is  shown  by 

defining  ^  as  in  (14b, c) ,  substituting  the  resulting  expression  for 
H 

(Oj^  +  K  into  (11a)  and  performing  the  average  over  the 

disturbance  ensemble. 

At  this  point,  we  note  that  solution  of  the  variational 
problem  of  (13)  and  (14)  requires  the  evaluation  of  the  ensemble 
average  of  p.  In  general  "p  satisfies  an  infinite  dimensional 
system  of  ordinary  differential  equations.  Thus,  for  practical 
use  this  system  has  to  be  closed  at  some  finite  stage.  For  this 
purpose,  general  results  are  given  by  Kistner.  There  appear 

to  be  two  situations  in  which  the  resulting  moment  equations  are 
tractable.  In  the  first  case,  the  Lie  algebra  generated  by 
(m  -6<)  and  v  is  solvable the  greatest  simplification  being 
obtained  for  an  Abelian  Lie  algebra.  The  second  case  is  that  in 
which  v(t)  is  white.  Since  general  conditions  on  the  Lie  algebra 
generated  by  (u  -6(c)  and  v  cannot  be  guaranteed,  the  strategy  we 
adopt  here  consists  in  replacing  the  actual  statistics  of  v(t)  by 
an  equivalent  (in  a  sense  to  be  discussed  below)  white  noise  model. 
Then  (13)  together  with  a  single,  closed  equation  for  ^  yield  a 


A _ . 


variational  problem  whose  solution  reduces  to  a  special  form  of 
the  well-known  state-dependent  noise  formulation. 

The  crux  of  the  matter  is  the  appropriate  choice  of  an  approx¬ 
imate  probability  model  for  v  (t) .  Before  pursuing  this  topic, 
we  first  state  some  preliminary  results  concerning  system  (11c, d) 
and  the  solution  of  (14b) . 

Theorem  1 

Suppose  that  v(t)  is  a  stationary  zero  mean  random  matrix 
process.  Define  an  increment  in  the  nonstationary  process  W  by: 


W(t^,t2)  k  W(t2)-W(tj^) 


di  V  (x) 


t2  ^  4 


(15) 


and  assume  that  W(tj^,t2)  is  almost  everywhere  continuous  in  tj^,t2 
and  bounded  for  all  finite  t2-tj^,  and  that  increments  of  the  form 
(15)  possess  joint  moments  of  all  orders.  Further,  suppose  that 
<(t)  is  bounded  and  continuous.  Then  with  p  and  3  as  defined  by 
(10)  and  the  foregoing  definitions: 


A.  The  transition  matrix,  <))(t,T),  for  system  (11c, d)  is  given  by: 


(t,T)  =  2  (J)  (t,T)  ;  t  >  T 

k=0 


a. 


where 


(t,T)  =  exp[u(t-x)  +  W(T,t)] 


b. 


>  (16) 


0.  (t,x)  =  (-l)*"/^ 

ilo) 


^k-1 


dx,/ 


dx  _ . . . / 

2  X 


dx, 


X  (<tiQ(t,Xi)  6k  (Xi)  (t»^  (Xi,T2)  (T2)  .  .  .4>o(Tk_i '■'^k^  ^^k^  ^  ^ 


J 


and  where  the  integrations  extend  over  the  left  semi-closed 
intervals . 
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A 


B.  ())(t,T)  is  almost  everywhere  continuous  and  its  first  and 
second  moments  are  continuous  and  differentiable  in  both 
arguments. 

C.  Eqs. (14b, c)  possess  the  unique,  positive  semi-definite 
solution: 


P  (t) 


t, 

/  ^  dT  tjj  (t , t) 
t 


t  e[tQ,tj^] 


(17) 


where,  for  t  2  t,  6s  x-t: 


<(»(t,x)  =  (t)*^(t+6,t)iJ;(t+6,T)<j)(t+6,t) 
i(j(t,x)  =  +  k^(t)R2<(t) 


(18) 


D.  tp(t,T)  4  E[i);(t,T)];  t  e  is  continuous  and  differen¬ 

tiable  in  t  and  x. 

The  proof  is  contained  in  Appendix  1.  At  this  point,  we 

remarJc  on  the  case  in  which  W(t)  is  a  matrix  Wiener-Levy  process. 

Expressions  (16)  through  (18)  retain  validity  but  the  resulting 

forms  of  averaged  quantities  are  not  implied  by  the  Ito 
(32) 

equation: 


dC(t)  =  (y  -  8K(t))C(t)dt  +  dv(t)C(t)  +  dw(t) 


(19) 


but,  instead  are  consistent  with  the  Ito  equations  with  the 
Stratonovit Jh  correction: 
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d^(t)  =  (y  -gK  +  I5I)  ^  (t)dt  +  C(t)  dv(t)  dw(t) 

I  t  LIM  i  E[W^ (t  t+A) ] 

A+0  ' 


y  (20) 


j 


In  other  words,  when  evaluating  response  averages  we  consider 
white  v(t)  as  the  limit  of  a  band-limited  process  as  the  band¬ 
width  approaches  infinity.  This  result  is  desirable  from  a 
practical  point  of  view,  since  parameter  deviations  actually 
encountered  are  likely  to  be  piecewise  continuous  (and  may  often 
be  random  variables  constant  in  time) .  Qualitative  features  of 
second  moment  response  which  obtain  for  v(t)  continuous  should 
be  preserved  when  v(t)  is  replaced  by  some  equivalent  white 
noise  model.  With  v(t)  piecewise  continuous  parameter  uncer¬ 
tainties  of  the  form  (12)  ,  it  is  easily  seen  from  (16)  that  the 
second  moment  response  of  the  uncontrolled  system  is  always 
asymptotically  stable.  While  this  result  is  duplicated  when  (20) 
is  employed,  (19)  implies  that  second  moment  stability  is 
dependent  upon  the  magnitude  of  the  white  noise  intensity. 

Thus,  in  the  white  noise  case,  use  of  (16)  through  (18)  or, 
equivalently,  (20)  is  indicated. 

2 . 2  Statistical  Modelling  of  Uncertainties 

As  Eqs. (17)  and  (18)  indicate,  a  complete  specification  of 
the  statistical  structure  of  open-loop  frequency  deviations  per¬ 
mits  explicit  determination  of  p,  so  that  the  optimization  scheme 
is  reduced  to  a  problem  in  the  calculus  of  variations.  However, 
complete  data  on  parcimeter  statistics  never  exists.  In  part  this 
arises  from  the  impossibility  of  devising  a  scheme  of  empirical 
inference  sufficiently  comprehensive  to  provide  valid  statistical 
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estimates  of  all  the  characteristic  functions  of  frequency  uncer^ 
tainties.  Further,  there  is  the  practical  difficulty  associated 
with  the  number  of  statistical  parameters  that  must  enter  into 
the  calculation  of  J.  The  number  of  statistical  measures 
(covariances  and  perhaps  higher  order  moments)  of  parameter 
uncertainty  must  perforce  increase  rapidly  with  the  order  of  the 
system.  For  a  very  large-order  system,  such  as  we  consider  here 
the  sheer  mass  of  uncertainty  measures  outdistances  our  ability 
to  enumerate  them  all. 

This  point  is  abundantly  illustrated  by  the  results  of 
Reference  (24)  wherein  the  stochastic  optimal  control  problem 
is  solved  for  a  general  multi-variable  discrete-time  system  with 
white  parameter  uncertainties.  The  resulting  Riccati-like  equa¬ 
tion  requires,  as  elementary  data,  the  covariances  of  several 

fully  populated  random  matrices.  For  each  such  matrix  this 

2  2 

entails  specification  of  %N  (N  +1)  scalar  covariances  (N  being 
the  system  dimension) .  Clearly,  for  N  large  the  difficulty  of 
interpreting  the  meaning  and  design  significance  of  all  these 
parameters  is  insuperable. 

Thus,  in  practice,  it  is  necessary  to  synthesize  insensitive 
controllers  given  incomplete  data  on  parameter  statistics.  As 
in  spectral  analysis  and  related  fields,  a  more  or  less  compre¬ 
hensive  probability  model  must  be  reconstructed  from  severely 
limited  data  in  a  manner  which  is  consistent  with  the  data  at 
hand  and  maximally  unpresumptive  with  regard  to  unavailable  data 
The  successful  principle  enunciated  by  Jaynes has  immedi¬ 
ate  application  here:  the  desired  probability  assignment  is  the 
one  which,  under  the  constraints  imposed  by  available  data,  can 
be  realized  in  the  maximum  number  of  ways  or,  equivalently, 
maximizes  the  entropy  of  the  underlying  processes. 


1 

A..  ....  . .  -  .  -  ■  ■ 


With  the  above  scheme  for  reconstructing  parameter  statis¬ 
tics,  it  is  clear  that  the  resulting  form  of  the  optimization 
problem  (i.e.,  the  form  of  p)  depends  critically  upon  the  nature 
of  the  available  data.  In  practice  it  may  be  advantageous  to 
acknowledge  as  accessible  a  set  of  data  which  is  even  more 
restricted  than  the  data  actually  available  but  which,  under  the 
maximum  entropy  principle,  renders  p  tractable.  In  the  follow¬ 
ing  we  attempt  to  discern  a  data  set  which  (1)  is  significant  to 
directly  observable  attributes  of  system  response,  (2)  consti¬ 
tutes  a  minimxam  set  for  the  purpose  of  maximum  entropy  probability 
assignments,  and  (3)  simplifies  calculation  of  the  parameter 
ensemble  averaged  performance  index.  ' 

First,  since  the  characteristics  of  modal  frequency  uncertain¬ 
ties  are  inherently  associated  with  the  uncontrolled  system,  let 
us  suppose  that  the  available  data  is  collected  from  measurements 
performed  on  realizations  of  the  uncontrolled  system  drawn  from 
the  parameter  ensemble.  Such  measurements  ultimately  entail 
direct  observation  of  the  attributes  of  system  response  to 
specified  disturbances.  For  example,  modal  frequencies  are 
never  directly  measurable  since  they  are  derived  quantities  which 
presuppose  a  dynamical  model.  Thus,  we  must  imagine  that  the 
available  data  consists  of  low-order  statistics  of  the  uncontrol¬ 
led  system  response  to  the  disturbance  noise  modelled  above. 

Since  practical  exigencies  would  preclude  estimation  of  higher 
order  moments  we  may  limit  consideration  to  first  and  second-order 
response  moments.  Clearly  the  expected  value  of  the  co-state 
matrix  must  also  be  considered  subject  to  direct  observation, 
particularly  as  its  value  for  the  uncontrolled  system  establishes 
the  fiducial  level  of  quadratic  cost. 

With  regard  to  these  quantities,  the  following  results  are 
easily  obtained. 

I 

j 

I 
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Theorem  2 


Consider  the  uncontrolled  system; 


^  =  (u  +  v(t))^+w(t) 


«"=o>  ■  io 


with  fixed  and  u ,  v  and  w  as  defined  previously.  Define 
increments  of  the  nonstationary  process  6^(0,t): 


=  ilm(M^)  6^  (ti  ,t.,) 


k'  ”k"'l'"'2' 


6^(t^,t2)  ^  6^(0, t2)  -6^(0, tj^)  ^  dT  6.  (t) 

t  ^ 


k  “  l/.«./2n  f  t2  -  tj^ 


with  the  log-characteristic  functions: 


^k^^'^1'^2^  =  E[exp[iu6^ (t^ ,t,) ] ] 


k'  l'^2' 


and  define  the  "modal  decorrelation  times",  Tj^; 


Tj^  £  (Ij^lm(uj^))  ^  di  |e 


ffc  dm  (Mj,)  ;  0  ,  T )  2 


)c  " 


where  the  I,  are  the  associated  reciprocal  time  constants. 


A 


Further,  assuming  the  Scime  properties  for  the  as 

in  Theorem  1  then  for  all  k;  t  e  [t^/tj^]  : 

2  2Reu-(t-t  ).  ^k  2 


A.  |E[5^(t)  ]  r  =  e“‘'-"k 


O'  e 


(24) 


«kk  -  ''kk4‘  k  e^R^'^kl'^-'o’ 

o 


t  (ii  +u*)  (t-x)  ( Imjjj^;  T  ,  t ) +r  •  (Imu  ■  ;  T  ,  t ) 

5kj  -  "kj^  "3’‘  '  ’ 

Mj  °  _  _ 


(25) 


(26) 


where  Q  denotes  the  second  moment  matrix  of 


Pkk(^) 


Pkj(t) 

(k?^j) 


(ai)kkP  dx 

^  ^  rS  /I  (y*+U.)  (x-t)  (ImUj  ;t,T)+rj^(ImMj^;t,x) 

(a,),  ./Idx  e  "^k  3  e  ■’ 

i  K]  t 


y  (27) 


LIM  Ip  I  P  |(Oi)kj|(T^T  ) 


(28) 


B. 


rkdroUk?  to'^^  I  2 

e  <1 


Tk'^  .  0 


a. 


b.  J 


(29) 


where  the  equalities  hold  if  5j^(t)  is  zero  almost  everywhere  for 
all  t. 
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Proof 

Result  (24)  follows  directly  from  averaging 


f.  (t) 


4>  (t,t^)C^ 
'  o ^ 


t 

/  4)q  ( t ,  T )  dw  (  T  ) 


with  if  (t,T)  given  by  (16b)  and  using  definitions  (21)  and  (22). 

^  r  H 

Similarly,  a  straight-forward  evaluation  o  f  E  ICC"]  yields  (25) , 
while  (26)  follows  by  use  of  the  Schwarz  inequality.  In  view  of 
(17)  and  (18)  with  <=0: 


(t)  =  /I  di  »t) 


Averaging  of  this  expression  leads  to  (27)  and  (28). 

Equation  (29a)  expresses  an  elementary  property  of  the  char¬ 
acteristic  function^^^^  while  (29b)  holds  by  definition  (23). 

Note  that  (29a)  ensures  second  moment  stability  and  the  existence 
of  constant  steady-state  values  of  ^  and  Q.  Now  if  S]^{t)  is  zero 
almost  everywhere  for  all  t,  then  also  vanishes  almost 

surely  for  all  Then  vanishes  and  the  equality  of  (29a) 

holds,  whence  the  integral  of  (23)  is  unbounded  and  the  equality 
of  (29b)  follows.  □ 
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As  (24)  shows,  the  mean  response  provides  the  magnitudes  of 
all  characteristic  functions.  However,  this  is  still  an  embarass- 
ment  of  riches,  and  we  must  seek  a  still  more  restricted  data  set 
v/hich  is  largely  independent  of  the  detailed  character  of 

In  this  connection,  (24)  and  part  B  of  the  theorem  reveal  an 
important  qualitative  effect.  Aside  from  small  natural  dissipa¬ 
tion,  the  energy  associated  with  the  mean  response  of  the  k^ 
mode  is  proportional  to  | e  |  which  tends  to  zero  as  t  tends  to 
infinity.  The  gross  effect  of  the  frequency  uncertainty  is  to 
introduce  a  spurious  damping  (the  "decorrelation  damping")  into 
the  mean  response  due  to  progressive  decorrelation  among  indivi¬ 
dual  ensemble  members.  Of  course,  as  (25a)  implies,  the  modal 
energy  thus  lost  to  the  mean  response  serves  to  augment  the 
covariance.  Thus,  the  mean  response  energy  gives  a  measure  of 
the  system  information  retained  at  any  time  subsequent  to  the 
application  of  a  known  disturbance. 

Now  the  ratio  of  the  mean  response  energy  to  that  predicted 
by  a  deterministic  model  is  the  magnitude  squared  of  the  char¬ 
acteristic  function  appearing  on  the  right  of  (24) .  This  ratio 
is  always  less  than  unity  and  integrable  so  that  any  modal 
frequency  uncertainty  renders  Tj^  finite.  Thus  the  decorrelation 
times  defined  by  (23)  are  the  time  scales  over  which  a  determi¬ 
nistic  model  retains  validity  given  known  initial  disturbances, 
and  appear  to  give  natural  measures  of  uncertainty. 

To  relate  the  Tj^  more  directly  to  frequency  statistics, 
suppose,  for  example,  that  is  known  to  be  a  zero-mean  Gaussian 
random  variable  with  standard  deviation  a.  This  corresponds  to  a 
relative  uncertainty  in  Im  with  variance  a  .  Using  (21) 

through  (23)  : 
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.  u 


or, 


and,  in  general  we  may  estimate  l/Tj^Im(yj^)  as  the  standard  devia¬ 
tion  of  the  frequency  deviation  relative  to  its  mean  value. 

Frequency  uncertainties  may  be  considered  "small"  when  the  1. 

iC 

are  small  or  when  the  Tj^  encompass  many  periods  of  natural  vibra¬ 
tion.  For  the  sufficiently  small,  the  modulation  introduced 
by  the  characteristic  functions,  exp  ( r,  (Im  (IT,  )  ;  t  ,  t)  )  ,  will  be 

JC  ic  o 

overwhelmed  by  the  attenuation  due  to  natural  damping.  This 
occurs  when  the  Tj^  are  much  larger  than  the  damping  time  scales, 
i. e.  : 

Tj^  >>  (Reilj^)"^  ;  k  =  l,...,2n  (30) 


In  this  regime,  as  can  be  seen  from  (24),  (25)  and  (27);  the 
influence  of  frequency  uncertainties  is  negligible  and  a  deter¬ 
ministic  model  of  the  plant  may  be  used. 

In  the  more  interesting  case  in  which  (30)  is  not  satisfied, 
various  qualitative  features  of  Q  and  p  may  be  deduced  from  (25) 
and  (27). 

For  this  discussion  and  henceforth,  we  suppose  that  the 
state  vector  5  is  so  arranged  that 


0)^  <  0.2  ^ 


S  0) 


n 
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and  define,  for  convenience,  the  integer- valued  function 
N (k) E [1 , 2 , . . . ,n]  such  that 


Further,  we  impose  various  restrictions  likely  to  be  satisfied 
in  practice: 


a. 


b. 

c. 


^2k”^2k-2  '  ‘^^^k'^k-1  with  increasing  k 


0  <  ^  n  ;  Vk 


Tj^llmUj^l  monotone  decreasing  with  increasing  k, 


i.e. 


Ik  >  Ik-1  ' 


Condition  (a)  is  satisfied  in  the  usual  case  wherein  modal 
frequency  separation  (considered  as  a  function  of  frequency)  is 
bounded  by  some  finite  power  of  the  nominal  modal  frequency. 
Condition  (b)  postulates  an  upper  bound  on  the  modal  damping 
coefficients.  (c)  states,  in  essence,  that  the  uncertainty  of 
open-loop  frequencies  relative  to  nominal  values  increase  for 
the  higher  order  modes.  Actually,  the  degradation  of  accuracy 
with  increasing  mode  number  need  not  be  monotone,  but  condition 
(c)  nevertnelesb  reflects  the  overall  trend. 
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Consider  now  the  steady  state  behavior  (t  +  0  and  tj^f°°  in  (25) 
and  (27) ,  respectively)  of  the  covariance  and  expected  cost 
matrices. 


Corollary  1 

Consider  the  model  order,  n,  arbitrarily  large.  Assume  con¬ 
ditions  (a)  through  (c)  above  and  denote  by  and  the 

results  obtained  under  a  deterministic  model  (r=0).  Then  for  the 
uncontrolled  system  in  the  steady  state: 

A.  With  the  smallest  integer  such  that: 


‘^N(k^)  - 


+^^k  -2  I 

c  c 


^  "= 


then  for  k>  k,j=k+2: 
c 


iQki 


'Z  Q.^Q .  . 

kk  j  j 


15k?' 


VqIWqIC) 

'^kk 


v/3— 7- 
^^kk^j  j 


I  Pkj 


vCTmTiST 

kk  33 


(31) 


a. 


^  (32) 


b. 


J 


B.  Similarly,  there  exists  a  smallest  integer,  k^j^  such  that 
(32)  holds  for  k  >  k  ^  and  j  e (k+2 , . . . ,k+2L) . 

CXi 

C.  Given  e  ,  >  0,  there  exists  a  k  >  k  sufficiently  large 

P  U  u  C 

that: 
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A 


Q 


j  j 


< 


£ 


Q 


kk^j  j 


< 


e 


P 


for  k  >  k^  and  all  j  k. 

Proof 

We  need  consider  only  Q  since  the  results  for  p  follow 
analogously.  In  the  steady-state: 


^kj 

k?<j 


^kk 


V.  ./  dx  e'^k  '  e  • 

o 


V  /”  dx  e^^®^k" 
^kk  o 


where  the  above  integrals  exist  by  condition  (b) . 

Clearly  the  diagonal  elements  of  Q  are  unaffected  by  frequency 
uncertainties  and  Qj^j^  and  identical.  From  (26)  ; 


|Q 


kj 


kj 


(Vj) 


while,  by  direct  calculation: 

Iv.  ■ 


loi°> I  . 


IKl 


IV^ji 


(33) 
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A 


SO  that : 


I  vT, +y  *  1 

=  k  3 ' _ 


^^k^k±2%‘^m±l 


where  k=2m  or  2m-l.  In  view  of  (31)  and  conditions  (a)  and  (c) , 
(32)  follows  directly.  Similarly,  by  virtue  of  (a)  and  (c) ,  we 
may  determine  a  smallest  integer,  k^2  such  that 


w. 


-A 

c2  c2 


"c2 


and  use  (26)  and  direct  calculation  of  to  show  (32)  for 

k  >  k^2'  3  e(k+2,  k+4).  Repetition  of  this  argument  proves  the 
assertion  of  part  B. 

Finally  part  C  follows  by  use  of  (26),  (28)  and  condition 
(c)  .  □ 

Note  that  parts  A  and  B  did  not  give  results  for  k  even  and 
j  odd  or  vice-versa.  However  in  this  case; 


l^^kj'  = 


(k)'^‘^N(j) 


=  o(— 

oU - i=_ 

''‘^N(k)'^‘^N(j 


/  dT  je 
o 


rj^(Imuj^;0,T)+rj  (ImUj;0,T) 


( 39 ) 

where  the  first  line  follows  from  Reimann's  lemma  ,  and  the  sec 
ond  by  the  Schwarz  inequality.  Thus,  from  (26)  we  may  estimate: 
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'kj 


-0/  — ) 


so  that  for  k  and  j  not  both  odd  or  even,  iQ.  . |  is  negligible  in 

j 

any  case. 

The  Corollary  shows  that  the  main  effect  of  frequency  uncer¬ 
tainties  is  to  suppress  cross-correlation  among  the  open-loop 
modes.  Indeed,  as  part  C  reveals,  the  portions  of  Q  and  0  cor¬ 
responding  to  sufficiently  high  order  modes  tend  to  become 
diagonalized  under  the  decorrelating  effect  of  uncertainties. 

These  results  allow  a  division  of  the  open-loop  modes  into 
various  qualitative  regimes.  An  important  line  of  demarcation 
is  provided  by  oj  as  defined  by  (31).  We  shall  term  the 
quantity  the  coherence  limit  in  frequency  since  it  locates 
the  onset  of  reduced  inter-modal  correlation. 

Suppose  k  and  j  exist  such  that  “'N  { j )  ‘^c*  Then  as 

one  sees  from  (33) : 


^kj 


so  that  in  this  quasi-deterministic  range  frequency  uncertainties 
have  little  influence. 

On  the  other  hand,  for  modes  much  above  the  coherence  limit, 
frequency  uncertainties  tend  to  obliterate  modal  cross-correlation. 
In  this  regime,  as  Corollary  1  shows,  there  exist  k,  ^suffi¬ 
ciently  large  that  the  correlation  coefficient  Qj^  j/ VQj^j^Q  j  j  is  as 
small  as  desired.  We  may  say  that: 
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OJ 

c 


defines  an  incoherent  range  in  which  the  open-loop  modes  are 
uncorrelated,  Q  approximately  diagonal  and  (from  (25a))  inde¬ 
pendent  of  parcimeter  statistics. 

Of  course,  the  same  general  behavior  can  be  deduced  for 
0.  In  particular,  the  sub-bloclc  of  p  corresponding  to  modes  in 
the  incoherent  range  is  approximately  diagonal  and  independent 
of  the  (Im(^j^);  o,t);  k=1...2n.  Thus  the  specific  form  of 
the  has  no  influence  on  that  portion  of  the  cost  contributed 
by  the  incoherent  range,  while  the  total  cost  is  primarily  depen 
dent  on  the  location  of  the  incoherence  limit.  We  conclude  that 
most  qualitative  features  of  open-loop  response  are  dictated  by 
the  magnitudes  of  the  Tj^  relative  to  the  other  time  scales  of 
the  problem. 

At  the  very  least,  we  must  require  that  any  approximating 
probability  model  of  the  k=1...2n  should  be  capable  of 

duplicating  the  general  behavior  described  above  -  in  particular 
it  should  preserve  the  time  scales  of  decorrelation  damping, 
provide  a  correct  estimate  of  the  coherence  limit  and  satisfy 
the  bound  given  by  (26)  for  the  cross-correlations  of  high  order 
modes.  This  is  possible  only  if  the  decorrelation  times  are 
admitted  as  fundamental  data. 

Thus,  henceforth,  we  propose  to  acknowledge  only  the 
Tj^  (k=1...2n)  as  the  "available"  data.  At  this  point  note  that 
in  practice: 
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so  that; 

62irj(t)  =  ~ 

Since  the  simplest  model  requires  use  of  the  (k=1...2n)  only, 
we  shall  acknowledge  this  relation  only  to  the  extent  that: 


T 


2m 


'^2m-l 


(34) 


With  this  choice,  it  remains  to  construct  a  full  probability 
model  which  presumes  as  little  as  possible  regarding  the  unavail¬ 
able  data.  In  other  words,  it  is  desired  to  determine  the  proba¬ 
bility  assignment  which  maximizes  the  entropy  of  the  processes 
5,  (t),  k=1...2n  subject  to  the  constraints  implied  by  (23)  and 
(34). 

In  general,  the  entropy  functional  to  be  maximized  should  be 
defined  on  the  joint  probability  distribution  functional  of  the 
processes  5j^(t),  k=l,...n  for  all  t  on  the  real  line.  To  avoid 
the  use  of  such  an  unwieldy  entity,  we  shall  proceed  heuristically 
and  define  the  maximum  entropy  probability  assignment  induced  by 
the  decorrelation  times  as  the  set  of  all  joint  statistics  of  all 
finite  sets  of  increments  (21)  obtained  from  the  solution  of  the 
following  problem  in  the  limit  as  T  and  N  tend  to  infinity: 

"Under  the  Restrictions; 


(a ) 


N  ^  ® 


( Imuj^ ;  0  ,mT/N) 


m=  1 

with:  T  ,  T 


2m 


=  T, 


'^2m-I  = 


=  1 , .  .  .  ,n 
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(b)  The  increments  >  0  are  stationary 

(c)  '5]^(t^,t2);  t^,t2  5  0  are  zero-mean  and  possess  finite,  non¬ 
zero  second  moments  for  |  t2-t |  >  0 

choose  the  joint  distribution,  F,  of  the  increments 
6j^(0,m^)  7k  =  l,*--,2n;m=l,***,N  (36) 

so  that  the  entropy : 

H  ^  -  ;dF  In  P 


IS  maximized." 

-where  P  denotes  the  probability  density  and  it  is  seen  that 
the  limit  of  the  sum  in  (35)  is  the  right  side  of  (23). 

Note  that  restrictions  (a),  (b)  and  (c)  differ  from  the  form 

usually  assumed  for  the  available  information  in  the  problem  of 
determining  the  maximum  entropy  probability  distribution.  ^ 

However,  in  the  limit  as  N  and  then  T  approach  infinity,  the 
solution  is  quite  simple  and  intuitively  plausible.  Its  deriva¬ 
tion  IS  contained  in  Appendix  2.  Here  we  summarize  the  conclusion 
as  f ol lows : 
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Theorem  3 

Assuming  that  the  processes 

<5],(0,t)  ^  6j^(T)dT  (37) 

k  1 f • a • f  2n 


possess  finite,  non-zero  variances  for  allte(0,<*>)  and  stationary 
increments,  the  maximum  entropy  probability  assignment  induced  by 
the  data: 


=  (Impj^Ij^)  ^  dT|e 


k  =  1  ,  a  a  a  ,n 


'^2m  ~  '^2m-l  '  ™  “  1».  .  •  ,n 


rj^(Impj^;0,T) 


>  (38) 


is  the  one  under  which  the  <5j^(o,t);  k=l, 
Wiener-Levy  processes  with  intensities  I 


aaa,2n  ats  independent 
k/^(k)  • 


E[6^(0,t) ] 


T 


‘^N(k) 


■2m 


=  I 


2m- 1 


m  —  l,aaa,n 


r  (39) 


With  the  as  available  data,  the  maximum  entropy  probability 
model  gives  independent  modal  frequency  uncertainties  -  a  property 
that  was  assumed  in  Section  2a  la  Furthermore,  this  statistical 
model  satisfies  the  assumptions  of  Theorem  1  and  those  results  may 
be  employed  directly  a 
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Note  that: 

r^(ImUj^;0,T)  =  -^^N(k)^kl^l 

Substitution  of  this  into  (24),  (25)  and  (27)  shows  that  all  the 
qualitative  features  noted  above  for  Q  and  p  are  preserved. 

In  this  connection,  it  is  convenient  to  compare  the  result 
obtained  for  ^  under  the  above  white  parameter  model  with  p  as 
computed  under  a  complete  statistical  model  having  the  same 
decorrelation  times.  As  a  practical  matter  we  must  limit  consid¬ 
eration  to  the  uncontrolled  system.  With  this  restriction, 

suppose  that,  in  reality,  the  open  loop  frequency  deviations,  6, , 

2  ^ 

are  normal  random  variables  with  variances  .  Then  by  appli¬ 
cation  of  (27),  the  steady  state  value  of  ^]^j(k?^j)  is  found  to  be 


7  -  77(G) 

^kj 


A  . 


.  ,-2  2  -2  2,  2 


(39)' 


and  the  diagonal  elements  of  "p  are  unaffected  by  frequency  uncer¬ 
tainties  and  need  not  be  considered  further.  Of  particular  con¬ 
cern  in  this  comparison  is  the  behavior  of  p  for  large  uncertain¬ 
ties  and/or  high  order  modes  (in  view  of  conditions  (a)  and  (c) ) . 
Thus,  consider  the  specific  case; 

v^k^k  +  >>  lUj+y^l 

Then  (39)a  yields  approximately: 
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-r(G) 


vT 

2 


Ikj 


(40)^ 


Now  the  decorrelation  times  may  be  computed  as: 


T 


k 


^°k‘^N(k) 


-1 


so  that  for  large  the  maximum  entropy  statistical  model 

yields : 


"kj 


^Ikj 

‘*'N(k)°k'^‘*^N(j)°j 


(41)^ 


Comparison  of  (40)a  and  (41)a  readily  shows: 


(42)^ 


where  equality  holds  only  if  o,  (i)„  .  =  a.  . 

■'  ■'  k  N(k)  j  N(3) 

Thus,  at  least  in  the  present  comparison,  the  maximum  entropy 
statistical  model  correctly  models  the  suppression  of  off-diago¬ 
nal  elements  of  p  due  to  frequency  uncertainties.  In  fact,  as 
(42) a  shows  the  model  of  Theorem  3  somewhat  underestimates  this 
diagonalizing  effect. 

Most  importantly,  the  are  modelled  as  white  noise  so 

that  posterior  learning  is  impossible  and  the  stochastic  control 
problem  is  nondual.  Although  physically  unrealistic,  the  white 
parameter  uncertainty  model  entails  an  extra  degree  of  caution 
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in  the  control  and  provides  a  worst  case  situation  from  the  point 
of  view  of  parameter  identification.  Indeed  the  model  may  be  used 
to  determine  performance  degradation  due  to  parameter  uncertainty 
and  to  assess  the  need  for  identification  and  adaptive  algorithms. 

2. 3  Determination  of  the  Optimal  Gain 

Let  us  now  resume  consideration  of  the  optimization  problem. 
Adopting  the  statistical  model  set  forth  in  Theorem  3,  we  obtain: 

Theorem  4 

With  the  maximum  entropy  probability  assignment  induced  by 
the  decorrelation  times  as  given  in  Theorem  3  and  K{t)  bounded 
and  continuous  in  t  e[t^,t^],  ^  as  defined  in  (14)  is  the  unique, 
positive  semi-definite  solution  of 

-p{t)  =  (t  -6K(t)  -l5J)^F(t) 

+  'p(t)(u  -6K(t)  -hi)  (40) 


+  J{p}  +  +  K“(t)R2K  (t) 


p(t^)  =  0 


where 


I  -  diag  [w^l2  /  /  •  >  • 


and  where,  for  any  square  matrix,  M: 


{M}  ^  '*^22  '  •  " '^2n,2n^ 
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Proof 

With  the  statistical  model  of  Theorem  3,  W(tj^,t2)  as  defined 
in  (21)  satisfies  the  restrictions  of  Theorem  1  and  the  results 
of  this  Theorem  may  be  used.  In  particular,  (18)  yields: 


)J;(t,T)  =  0^  (t+6  ,t)  (t+6  ,t)  4)  (t+6  ,t) 

0(1+6, t)  =  0^  (t  ,  t+6)  [CTj^+fc^(T)  R2<  (t)  ]  0  (T  ,  t+6) 


a. 


b. 


(43) 


where  6>o.  (16)  shows  that  0{t  +  6,t)  depends  upon  '5]^{tj^,t2) 

only  for  t^,t2  (t,  t  +  <5],  while  0(t,  t  +  5)  depends  upon 

‘^k^^l'^2^  only  for  tj^,t2  e  (t,  t  +  6].  Since  these  intervals 
are  disjoint  and  the  increments  of  <5j^(o,t)  are  independent, 
the  ensemble  average  of  (43a)  becomes: 

0(t,T)  =  E[0^(t+6,t)0(t+6,T)0(t+6,t) ]  (44) 


where 


0(t,T)  =  E[0(t,T) ] 


NOW  examine  (16).  Keeping  in  mind  that  (1)  the  0q's  appear¬ 
ing  in  the  integrals  of  (16c)  are  each  dependent  on  increments 
of  the  6,  (o,t)  over  mutually  disjoint  intervals,  (2)  that  partial 
sums  of  (16a)  are  almost  everywhere  convergent,  and  (3)  the 
6,  {o,t)  are  Gaussian  with  variances  (39),  it  is  seen  that  the  con- 

K  QO 

tribution  of  T 
2 

of  order  6  on  the  right  side  of  (44)  . 


(t  +  <5,  t)  to  0(t  +  <S ,  t)  produces  terms  of 


(44)  may  thus  be  written: 
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.p(t,T)  =  E[  (<!>Q(t+6,t)  -/J^'^^dTj^(t)Q(t+6,T^)BK(Tj^)(l>Q(Tj^,t))” 

X  ?  (t+5  ,  t)  (4)q  (t+6  ,t)  -/^^'*’‘^dTj^())^(t+6,Tj^)  Btc  (Tj^)  (})^(T^,t)  )  ] 

+  0(6^) 

Similarly,  using  tVe  expression  (16b)  for  the  appearing 

above,  we  have; 

=  E[A^^(t+6,T)  A]  +  0(6^) 

A  =  I  +  u6  +  W(t,t+6)  +  ^W^(t,t+6)  -/^^■^‘^dT^6<  (Xj^) 

After  expanding  out,  rearranging  and  dividing  by  6: 

-  ■^[4)(t+6  ,t)  -4)(t,T)  ] 

=  (M-6K(t)  +  ^E[W^(t,t+6)  ]  )^ljJ(t+6,T) 

+  ^(t+6,T)  (ir-6<(t)  +  ^E[W^(t,t+6)  ]) 

+  jE[W^(t,t+6)iiJ(t+6,T)W(t,t+6)  ]  +  0(6) 

Next  use  (21a)  and  Theorem  3  to  evaluate  the  above  averages, 

9  — 

then  pass  to  the  limit  64-0.  Recalling  that  ij;  (t,T)  exists  by 
Theorem  l.D,  we  obtain: 

-  ^^(t,!)  =  ()a-6K(t)  -hl)^  ■iF(t,T) 

+  'i(j(t,x)  Cu-Bk  (t)  ~hl)  +  l{^(t,T)}  (45a) 

with  I  given  as  in  (41).  Also,  (18b)  yields  directly: 

ip(T,t)  =  +  k^(t)R2k:(t)  (45b) 
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Finally,  integration  of  all  terms  in  (4  5)  over  t  e  [t,tj^]  and 
use  of  (17)  gives  (40). 

The  linearity  of  this  equation  guarantees  the  uniqueness  of 
the  solution,  and  the  positive-semidefiniteness  of  p(t)  noted  in 
Theorem  l.c  implies  the  seune  property  for  ^(t).  □ 

Under  the  maximum  entropy  statistical  model  we  thus  obtain 
a  modified  Lyapunov  equation  for  p  which  must  be  appended  to  the 
variational  problem  (13)  as  a  constraint.  Clearly  as  the  decor¬ 
relation  times  approach  infinity,  the  matrix  I  approaches  zero 
and  (40)  reduces  to  the  familiar  Lyapunov  equation  for  a  deter¬ 
ministic  plant. 

With  finite  decorrelation  times,  the  qualitative  structure 
of  (40)  should  be  noted.  First,  considering  only  the  diagonal 
elements  of  (40) ,  we  have 

-  4{p}  =  u”{p}  -{(6k)”?} 


+  {ply  -{p6k}  +  {oj^+k”r2k} 


Thus  the  terms  arising  from  frequency  uncertainties  do  not 
appear  explicitly  and  the  diagonal  elements  of  ?  depend  on  the 
decorrelation  times  only  through  their  dependence  upon  the  off- 
diagonal  elements. 

On  the  other  hand,  considering  only  off-diagonal  elements, 
(40)  takes  the  form; 


-p  =  (y-6K-^I)”p  +  p(y-6K-*iJ) 

+  +  k^R2'^ 
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The  off-diagonal  elements  of  P  are  directly  influenced  by  the 
uncertainty  terms,  but  only  through  the  matrix  (u"  -6k  -hi)-  It 
is  easily  shown  that  under  the  statistical  model  of  Theorem  3, 
the  mean  respoiise  (averaged  over  the  parameter  ensemble)  ,  Xi  is 
given  by 

'C  =  (p-gK-Js!)!’ 

Consequently  (u  -6k  -hi)  is  the  dynamic  matrix  of  the  "mean 
system"  and  the  term  -hi  represents  the  effect  of  decorrelation 
damping. 

Thus,  while  the  diagonal  elements  of  "p  are  not  directly 
affected  by  parameter  uncertainties,  the  off-diagonal  elements 
are  subject  to  the  decorrelation  damping  of  the  mean  system. 

These  features  lead  to  the  same  suppression  of  off-diagonal  ele¬ 
ments  of  Q  (and  Q)  as  noted  previously  for  the  uncontrolled 
system. 

The  problem  outlined  in  (13)  and  (14)  is  now  reduced  to  one 
in  the  calculus  of  variations.  With  (13)  and  (40)  we  may  pro¬ 
ceed  directly  to  obtain: 

Theorem  5 

Under  conditions  (40),  the  performance  index: 

—  ^1  — 

J  =  /  dt  tr[pv]  (46) 

^o 

is  minimized  for  all  v  ^  0  if  and  only  if: 

K  =  (47) 

where  ^  is  the  positive  semi-definite,  hermitian  solution  of 


which  may  be  recognized  as  the  equation  determining  the  second 
moment  matrix.  Equating  to  zero  the  first  variation  of  (50) 
with  respect  to  k  yields; 


(R2K  -6”p)Q  =  0 


For  this  to  hold  for  all  v,  (47)  must  follow.  Finally,  substi¬ 
tution  of  (47)  into  (40)  produces  (48). 


-p  =  (y-i5i')^p  +  p(y-i5i)  +  I{p} 


+  a^-pa2P 


P(t^)  =  0 


where 

02  =  (49 

Proof 

The  necessary  stationary  conditions  may  be  derived  by  intro¬ 
ducing  a  multiplier  matrix,  Q,  and  requiring  that  the  first 
variation  of 

^1  _  H  —  _  H _ 

H  =  /  dt  tr[pv  +  Q"[p  +  (u-BK-*5f)"p 

(50 

+  'p(u-6k->5J)  +  I{^}  +  Oj^  +  k^R2<]] 

with  respect  to  independent  variations  in  "p  and  <  vanish,  impos¬ 
ing  also  the  terminal  condition  (40b) .  Vanishing  of  the  first 
variation  with  respect  to  "p  gives; 


Q  =  (m-Bk-JsDQ  +  Q(y-6K-»iI)"  +  J{Q}  V 


Qltg)  =  0 
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To  prove  sufficiency  of  (47)  and  (48),  first  define: 


K 


o 


where  satisfies  (48).  In  general,  we  may  set 

K  =  K„  +  K 
O 

P  =  Pq  +  z 

With  this  substitution,  (40)  becomes: 

-z  =  (]I-^I- Sk)  ^2  +  z(y-55l-6<)  +  l{z}  +  k^R2k: 
z  (tj^)  =  0 


(52) 


Next,  define  q(t)  by: 

q  =  (u-m-6K)q  +  q(y->5l-Bic)  +  l{q} 

qCt^)  =  %  >  0 

where,  by  construction,  q  >  0  on  t  e  [tQ,t^]. 
(53)  : 


(53) 

Then,  with  (52)  and 


^  tr[zq]  =  -tr[<^R2Kq]  <  0 


V  t 


(54) 


where  the  inequality  is  implied  by  q  >  0,  R2  >  0.  From  the  ter¬ 
minal  condition  on  z: 

(tr[zq])  =  0 

^1 
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Combined  with  (54),  this  implies: 

tr[zq]  >  0  ;  t  e{t^,tj^] 
and  since  q  >  0: 

P  -pQ  s  0  ;  V  t  c[t^,tj^],  <9^0 
Therefore,  since  v  >  0: 

J(k)  >  J(<^) 

which  shows  that  J  is  minimized  with  (47)  and  (48) .  D 

So  far  the  requirement  that  K  be  real  has  not  been  mentioned. 
But  this  condition  is  inconsequential  in  view  of  the  following 
result: 

Corollary  2 

Let  p  be  any  hermitian  solution  of  the  terminal  value  problem 
(48).  Then  the  matrix,  P,  of  iP  expressed  in  the  modal  coordinate 
basis : 

P  ^  (55) 

is  real  and,  consequently, 

K  =  R2^  P  (56) 

is  real. 

The  elementary  but  laborious  proof  is  given  in  Appendix  3. 

As  a  consequence,  subsequent  theoretical  developments  may  be 
carried  out  within  the  eigen~basis  of  A,  employing  the  complex 
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equation  (48).  In  numerical  work,  however 
advantage  to  present  (48)  in  the  original 
Use  of  (55) ,  (48)  and  previous  definitions 

-P  =  (A->5.")'^P  +  P(A-l5/)  +  D(/,P1 
+  -P  B  r“^b'^P 

P(tj^)  =  0 


,  it  is  of  obvious 
modal  coordinate  basis, 
gives : 


(57) 


where,  for  real  P; 


D[/,P] 


block-diag 

k=l , . . . , n 


^^2k-l  ,2k-l'^“k^2k,  2k^ 


2w 


oj  0 

k 

0  1 


(58) 


This  expression  for  D[I,P]  serves  to  illustrate  the  relative 
simplicity  of  the  complex  form  given  by  (48) . 

2 . 4  Recapitulation  of  the  Proposed  Approach 

In  addressing  the  problem  of  optimal  control  of  uncertain 
structural  systems  in  the  preceding  sections,  we  have  developed 
an  approach  which  is  an  explicit  expression  of  a  general  design 
philosophy  which  differs  significantly  from  the  traditional  view. 
At  this  point  it  is  well  to  contrast  the  proposed  formulation 
with  conventional  application  of  modern  control  theory  and  out¬ 
line  the  ramifications  of  the  new  design  strategy. 

Such  a  comparison  is  depicted  in  general  terms  by  Fig.  1. 
Clearly,  the  traditional  approach  has  a  rich  and  diverse  back¬ 
ground,  yet  we  must  paint  it  with  broad  strokes  in  order  to 
reveal  its  essential  premises. 
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T  H  - 


Conventional  Alter  n a  t  e 

Philosophy  Ap  p  r  o  a  l  m 


Fig.  1.  Stochastic  optimal  control  under  limited  s\'st(3n  information. 
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Thus,  as  indicated  on  the  right  of  Fig.  1,  the  alternate 
approach  begins  with  a  severely  limited  set  of  parameter  data. 

In  practice,  we  are  limited  to  the  nominal  or  mean  values  of 
system  parameters  and  to  statistical  measures  of  the  variation 
about  mean  values.  There  is  some  latitude  of  choice  with  regard 
to  measures  of  parameter  variation  and  it  is  expeditious  to 
acknowledge  as  available  a  data  set  which  is  also  significant 
to  the  overall  fidelity  ^f  the  model. 

Having  identified  a  significant  and  limited  data  set  which 
we  choose  to  acknowledge  as  available,  a  complete  probability 
model  is  still  required  for  computation  of  the  control  law. 
Consistency  requires  that  we  induce  this  complete  probability 
model  uniquely  from  the  acknowledged  data.  To  do  this  we  follow 
the  doctrine  that  the  complete  statistical  assignment  be  consis¬ 
tent  with  acknowledged  data  and  maximally  unconstrained  other¬ 
wise,  i.e.,  we  appeal  to  a  maximum  entropy  principle.  This 
introduces  an  element  of  design  conservatism  since  the  statis¬ 
tical  "degrees  of  freedom"  of  the  maximum  entropy  model  are  far 
greater  in  number  than  under  the  actual  parameter  statistics. 

Thus  stability  properties  found  under  the  maximum  entropy  model 
are  very  likely  to  hold  over  the  actual  parameter  ensemble. 

At  this  stage  (with  specification  of  items  1  and  2  in  Fig.  1) 
we  possess  the  "system  model"  employed  by  our  alternate  approach. 
This  model  is  necessarily  of  high  dimension,  as  is  the  truth 
model,  but  implicitly  accounts  for  the  substantial  uncertainty 
in  high  order  modal  parameters  while  requiring  relatively  little 
elementary  data  as  input. 

With  the  complete  probability  model,  we  must  compute  the 
control  law  based  upon  a  performance  measure  defined  on  the 
entire  parameter  ensemble.  The  simplest  choice  appears  to  be 
the  average  of  a  quadratic  criterion  over  the  parameter  ensemble 


since  it  pro/o-des  us  with  a  straightforward  mechanism  whereby 
almost  sure  stability  may  be  guaranteed. 

Under  the  restrictions  of  full-state  feedback,  and  uncertain¬ 
ties  only  in  the  open-loop  frequencies,  the  preceding  sections 
have  given  specific  form  to  items  1  through  3  indicated  by  Fig.  1. 

To  recapitulate,  we  first  attempted  to  identify  a  data  set  which 
was  significant  to  modelling  fidelity  by  examining  the  phenome¬ 
nology  of  frequency  uncertainties  as  reflected  in  the  mean  response, 
the  covariance  and  the  expected  cost.  As  a  result,  we  identified 
the  modal  decorrelation  times  as  the  parameter  data  set  that  must 
be  acknowledged  as  available  if  the  open-loop  system  second  moment 
response  is  to  be  adequately  modelled.  Indeed,  various  qualita¬ 
tive  features  such  as  (1)  the  location  of  the  "coherence  limit," 

(2)  the  behavior  of  the  portions  of  the  expected  cost  matrix  and 
second  moment  matrix  pertaining  to  modes  in  the  "incoherent  range," 
and  (3)  the  gross  effect  of  "decorrelation  damping"  depend  explic¬ 
itly  upon  the  decorrelation  times.  Choosing  to  acknowledge  only 
these  quantities  as  available  data,  we  constructed  a  full  proba¬ 
bility  model  for  frequency  uncertainties  which  is  otherwise  maxi¬ 
mally  unconstrained.  The  resulting  statistical  model  is  an 
equivalent  white  noise  parameter  model  which  permits  formulation 
of  closed  equations  for  the  expected  cost  matrix  and  explicit 
solution  of  the  linear  regulator  problem. 

One  clear  advantage  to  the  present  formulation  is  the  rela¬ 
tively  small  number  of  parameter  statistical  measures  that  must 
be  provided  and  their  direct  relationship  to  system  characteris¬ 
tics.  Secondly,  the  white  parameter  model  yields  a  non-dual 
problem.  Thus,  while  accepting  frequency  uncertainties  at  their 
2  rr‘'or{  levels,  the  approach  is  consistent  with  the  more  general 
context  of  stochastic  optimal  control. 


In  addition,  certain  features  arising  from  the  specific  form 
of  (48)  are  essential  to  the  overall  rationale  and  will  be  eluci 
dated  in  the  sequel.  To  anticipate  these  results,  it  will  be 
shown  that  under  weak  restrictions,  (48)  possesses  a  unique  posi 
tive  semi-definite  solution  in  the  steady-state  case  for  all 
values  of  the  decorrelation  times.  This  precludes  the  existence 
of  an  uncertainty  threshold.  Also,  with  the  gain  (47),  second 
moment  stability  is  guaranteed  for  the  closed-loop  system.  Thus 
the  need  for  design  iteration  to  ensure  robustness  with  respect 
to  stability  is  largely  eliminated  within  the  present  approach. 

Finally  the  form  of  the  modified  Riccati  equation  has  impor¬ 
tant  implications  for  the  dimensionality  issue.  We  shall  find 
that  by  confronting  parameter  uncertainties  directly,  we  shall 
also  greatly  reduce  the  need  for  modal  truncation  with  its 
attendant  spillover  problems. 
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3.  LINEAR  REGULATION  UNDER  MODAL  FREQUENCY  UNCERTAINTIES: 
RAMIFICATIONS  OF  THE  STOCHASTIC  DESIGN  APPROACH 


3 . 1  Introduction 

In  the  remainder  of  this  report  various  ramifications  of  the 
design  approach  outlined  previously  are  investigated.  Specifi¬ 
cally,  we  resume  consideration  of  linear  full  state  feedback 
regulation  in  the  presence  of  open-loop  frequency  uncertainties 
and  develop  the  properties  of  solutions  to  the  terminal  value 
problem  (48) . 


First,  we  recapitulate  the  formulation  and  specify  the 
restrictions  under  which  the  work  will  proceed.  It  must  be 
noted  that  previous  results  retain  validity  if  =  0  in  (2)  for 
some  k,  i.e.,  if  rigid  body  modes  are  explicitly  recognized.  In 
the  following  the  specification: 


^2nx2n 
c  C 


(59) 


n  4  n^  +  n^ 

-  r  e 

will  be  assumed,  where  is  the  dynamic  matrix  of  n^  rigid  body 
degrees  of  freedom: 


1^  =  block-diag 


0  1 
0  0 


and  Ug 


encompasses 


the  elastic 


odes : 


(60) 


Ug  £  diag{(jj^  (i-rij^)  f  w.  {-i-n^) 


•■•“n  '^'V'  “n  '> 

e  e  e  e 


0  <  (jj,  <  w_,  <  ...  <0) 

12  n 


although  the  results  are  readily  generalized,  we  assume  for 
simplicity  that  all  the  are  distinct.  Further,  non-zero 
structural  dcimping  will  be  assumed  for  all  elastic  modes: 


>  0  ;  k=l , .  .  .  ,ng 


Consistently  with  (59),  the  matrix  of  inverse  decorrelation 
times  now  assumes  the  form: 

^[00] 

[o  leJ  .  (6 

le  =  diag[wj^l2,wj^l2r- •  •  l2^  12^  ] 

e  e  e  e 


with  the  as  defined  in  (38)  . 

Under  conditions  (59)  through  (63),  we  consider  the  termi¬ 
nal  value  problem  arising  from  Theorem  5; 


50 


-P  =  u“p  +  p  Mjji  +  +  aj^-pa2p 

p{t^)  =  0 


(6 


where 


-  ^  -  VT 


02  =  Br'^b”' 


>  0 


R2  >  0 


S  ^ 


e  J 


^  = 


rll 


rl!l 


L^rn  1 
r 


J 


(6 


{6 


(6 


and  where  e  is  as  defined  by  (lOe)  and  (6b)  and  B^  is 


as  defined  by  (10b) ,  (8)  and  (3) 
With  the  control  gain: 


K  =  R2^B^P 


(6 


the  second  moment  matrix  of  the  closed  loop  system  is  determined 
by : 

Q  =  (a  -hT)Q  +  Q(a  -hl)^  +  I iQ)  +  v 
Q{tQ)  =  >  0  ,  V  >  0 


where 


(6 


r 


A 


as  (51)  in  the  proof  of  Theorem  5  shows.  Here,  we  assume  that 
the  initial  state  is  non-zero  and  statistically  independent  of 
disturbance  and  parameter  noise. 

The  next  section  reviews  various  concepts  and  results  con¬ 
cerned  with  stochastic  stability.  Sufficient  conditions  for  the 
existence  and  uniqueness  of  steady  state  solutions  to  (64)  are 
established  in  Sections  3.3  through  3.5.  Much  of  this  material 
parallels  corresponding  developments  in  Refs.  (40)  and  (41) , 
except  that  several  restrictions  previously  stated  can  be 
removed  for  the  problem  considered  here.  In  addition.  Section 
3.5  achieves  assurance  of  various  stochastic  measures  of  stabil¬ 
ity.  For  the  steady  state  case.  Section  3.6  gives  certain  well 
known  numerical  procedures.  Finally,  the  section  concludes  with 
consideration  of  asymptotic  properties  for  high  uncertainty 
levels  and  high  order  modes  in  Sections  3.7  through  3.9. 

3. 2  Stochastic  Stability  and  Extended  Lyapunov  Equations 

Preliminary  results  concerning  the  stability  properties  of 
the  closed-loop  stochastic  system  (1)  must  be  given.  Any  bound¬ 
edness  or  convergence  property  used  in  deterministic  system 
theory  can  be  translated  for  stochastic  systems  in  different 
ways,  depending  on  the  type  of  stochastic  convergence  one  wishes 
to  consider. 

The  convergence  of  sample  solutions  of  a  stochastic  dynamic 

system  of  the  form  (1)  to  the  null  solution  are  most  frequently 

( 42 ) 

characterized  by  the  following  stability  definitions: 

1.  Almost  sure  exponential  stability:  The  sample  solutions 
converge  exponentially  to  zero  as  t  tends  to  infinity  with 
probability  one. 
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2. 


th 
2. _ 

and 

2n 

2 

k=l 


moment  exponential  stability:  With  p  a  positive  integer 
Pj^...P2j^  any  set  of  non- negative  integers  such  that 

Pk  =  P 


and  for  all  initial  states. 


the 


p‘'' 


moments : 


(t) . .  .X 


2n 


(t) 


P2n 


] 


converge  exponentially  to  zero. 

3.  mean  exponential  stability:  For  p  >  0  and  all  initial 

states , 

E[  I  |x(t)  I  |P] 

converges  exponentially,  where  |1...1|  denotes  some  vector 
norm. 

Clearly,  p^^  moment  exponential  stability  is  useful  as  a 
convergence  property  only  if  p  is  an  even  integer  and,  with  p 

even,  definitions  2  and  3  are  actually  equivalent.  Moreover, 

/N,  th  '^th 

for  p  >  p  >  0,  p  mean  exponential  stability  implies  p  mean 

exponential  stability  and  almost  sure  stability. 

Since  the  system  performance  is  characterized  by  a  quadratic 
functional,  it  suffices  for  our  purpose  to  consider  only  second 
moment  (or  equivalently  second  mean)  exponential  stability.  Obvi¬ 
ously,  for  second  moment  exponential  stability  it  is  necessary 
and  sufficient  that 


lim  Q(t)  =  0  (70) 

t- 1  +«> 
o 

for  all  Qq  >  0  and  v  =  0  in  (69) .  If  the  system  is  second  moment 
stable  then  Q(t)  with  v  >  0  has  a  unique  equilibrium  value 
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X 


(71) 


lim  Q(t) 
t- 1 

o 


which  from  (69)  is  the  positive  semi-definite  solution  of 


0  =  X(a  +  (a  -hI)X  +  /{X}  +  v 

V  s  0 


y  (72) 


Particularly  germane  to  the  existence  of  steady-state  solu¬ 
tions  of  the  stochastic  Riccati  equation  (64)  is  the  character 
of  solutions  to  the  adjoint  of  (72) : 


0  =  A(a  -hi)  +  (a  -hl)^h  +  /{M  +  S 
S  >  0 


(73) 


(43) 

Conditions  previously  given  for  the  existence  and 

uniqueness  of  positive  semi-definite  solutions  to  stochastic 
Lyapunov  equations  of  the  form  (72)  or  (73)  may  be  stated  in 
terms  of  equivalent  coefficient  matrices  defined  as  follows: 
Since  Q  is  hermitian  we  may  set  up  the  vector 


q”  =  /2  Q22' 


(74) 


so  that  only  the  upper  triangular  portion  of  Q  is  used  and 
Q  e  ^n(2n+l)^  Vector  v  is  similarly  defined  on  v  in  (69). 
Then  (69)  may  be  written: 
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(75) 


A  A  A 

Q  =  AQ  +  V 


where  A  results  from  representation  (74)  and  is  termed  the 
"equivalent  coefficient  matrix  of  Q  in  [Q  (a-^J)  (a-^sl )  Q+7  (Q)]  •" 
To  display  its  definition  explicitly,  we  denote  A  in  (75)  by 


AQ[Q(a  -hl)^  +  (a  -hDQ  +  KQ}] 


Proceeding  similarly  for  A,  Equation  (73)  may  be  written; 
0  =  A^[A(a  -HI)  +  (a  -HD^ h  +  J{A}]A  +  S 


Clearly,  from  the  above  definitions; 


(76) 


and  a  necessary  and  sufficient  condition  for  second  moment  expo¬ 
nential  stability  is  that  Aq  (or  equivalently  be  asymptoti¬ 

cally  stable. 
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The  relations  between  the  stability  of  and  the  character 
of  solutions  of  (72)  and  (73)  may  be  illustrated  by  the  follow¬ 
ing  two  results.  First  we  have: 

Theorem  6 

Consider  (73)  with  S  ^  0  (S  s  0)  .  If 
A^[A(a  -hi)  +  {a  -hl)^A  +  I(A}] 


is  asymptotically  stable  then  there  exists  a  unique  positive 
semi-definite  (negative  semi-definite,  resp. )  solution  to  (73) . 

Proof 

Consider  the  case  S  >  0  first.  For  any  e  C^"^ ,  define 
Q(t)  by  (69)  with  v  =  0,  =  0,  and 

Q  ^  ^ 

O  -  ^0^0 


A 

SO  that  Q  is  determined  by 


Q  =  Aq  Q 
A  A 

Q(0)  =  Qq 


(77) 


If  is  asymptotically  stable,  then  so  is  A^  by  virtue  of  (76) . 
Consequently  (77)  possesses  the  unique  solution; 


Q(t)  =  e  ^ 
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with  Q(t)  positive  semi-definite  on  [0,  ^) . 


Similarly,  with  asymptotically  stable, 

A  A 

0  =  A, A  +  S 


has  the  unique  solution: 


A 

A 


oo 

/  dt 


With  these  results  we  have: 


H  oo  ^ 

A  =  tr{Q„A}  =  q“a  =  /  dt  Q“(t)S 

O  O  O  O  Q 


Furthermore,  the  inequalities: 


Q^(t)S  =  tr{Q(t)S}  >  I  |Q(t)  I  I  2  ?  0  (78) 

hold  for  t  e  [0,  “) ,  where  A{M}  denotes  an  eigenvalue  of  matrix 
M  and  1 1... 1 I2  denotes  the  spectral  norm.  Hence: 

A  °  ^o  ^  0 


and  the  assertion  for  S  >  0  is  shown. 
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For  S  ^  0,  we  may  replace  (78)  by 

Q^(t)S  =  tr{Q(t)S}  <  2n  S )  1  1  Q  ( t )  1  1  ^  ^  0 


so  that: 

«o  ^  «o  ‘  0  ■  'O  “ 

which  completes  the  proof.  D 

To  a  considerable  extent,  the  converse  of  the  above  may  be 
shown: 

Theorem  7 

Suppose  that  (a-^i)  is  asymptotically  stable,  S  >  0  and 

_  k 

(a-*sI,S^)  is  completely  reconstructible.  If  (73)  possesses  a 
positive  definite  solution,  then 

A^[A(a  -hi)  +  (a  -hl)^h  +  KM] 


is  asymptotically  stable. 

Proof 

By  Lemma  3.1  of  Reference  (44),  complete  reconstructibility 
—  h 

of  icL-hl  ,S^)  implies  that 


r  dt 


s 


^(a-hl )t 
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IS  positive  definite.  With  this  condition  and  the  remaining 
assumed  conditions,  the  stated  conclusion  follows  by  application 
of  Theorem  5-3  of  Ref.  (41). 


3 . 3  The  Time-Dependent  Equation  -  Convergence  to  the 
Steady-State 

Here  we  investigate  certain  properties  of  (64)  with  a  view 
to  establ ishi -ig  conditions  under  which  the  solution  approaches 
a  steady  state  value  as  t^^  tends  to  infinity. 

The  first  preliminary  result  allows  the  problem  to  be 
reduced  to  a  consideration  of  the  completely  reconstructible 
portion  of  the  system. 


Lemma  1 


A. 


Under  the  conditions  specified  in 
-  % 

(u  ,0  ')  is  detectable.  Then  the 
m  l 

rearranged  to  assume  the  form: 

r«  1 


so  that  E,  contains  the  rigid  body 
the  unreconstructible  subspace  of 
matrices  accordingly)  : 


Section  3.1,  suppose 
state  vector  may  be 


modes,  comprises 

and  (partitioning 


(79 


all 


A  A  A 

=  u  -hi 
m 


>  (8 


where  is  completely  reconstructible  and  I  is  a 

diagonal  matrix  comprising  those  entries  of  /  correspond¬ 
ing  to  the  reconstructible  portion  of  the  system. 
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B.  With  all  matrices  partitioned  in  a  manner  consistent  with 
(79),  and  with  conditions  (65)  and  (66),  the  solution  of 
(64)  exists,  is  positive  semi-definite  on  (-“,  t^]  and 
assumes  the  form: 


P  = 


(811 


where  p  is  the  solution  to 


d  -  -H- 

p  =  U^P  t  P 


p  ( t  ^ )  =  0 


7{p}  +  -pa2P 


:82) 


where : 

A 

A  / 

8 


m:.) 


u 

02  =  6R2  S 


(83) 


The  proof  is  contained  in  Appendix  4. 

Now,  upon  consideration  of  (82),  Theorem  6-5  of  Reference 
(41)  yields: 

Lemma  2 

PI  ^  Jj 

Consider  (82)  and  (83)  with  (u  ,0,  )  completely  recon- 

m  i 

In  addition  to  the  conditions  of  Lemma  1,  suppose 


structible. 


has  a  unique  positive  semi-definite  solution,  and  that 


i{Q)) 


is  asymptotically  stable.  Then  p(t)  >  0  for  all  t  <  t^^  and: 

A 

lira  p (t)  =  (85) 

t^  f  00 

Further,  with  the  reduction  afforded  by  Lemma  1,  conditions 
for  second  moment  stability  can  be  simplified  as  follows: 

Theorem  8 

Under  the  assumptions  of  Lemma  2,  and  denoting  the  steady- 
state  solution  of  (64)  by  p;  then  both  and: 

AQ[(Um  -a2^)Q  +  Q(Uj„ 

are  asymptotically  stable. 
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The  lengthy  but  elementary  proof  is  relegated  to  Appendix  5. 

From  the  above  result  and  Lemma  2,  the  situation  with  respect 
to  convergence  of  c  to  a  steady  state  solution  and  overall  second 
moment  stability  under  the  resulting  steady  state  control  law  is 
seen  to  hinge  on  the  existence  and  uniqueness  of  a  positive  semi- 
definite  solution  of  (84)  for  which  Aq  [  (iJ^-a2A)  Q+Q  (Pj^~0  2A)  {Q}  ] 

is  asymptotically  stable.  This  matter  is  addressed  in  the  next 
section . 


3 . 4  The  Time- Independent  Equation:  Existence  and 
Uniqueness  of  Solutions 

Here  we  determine  sufficient  conditions  for  the  validity  of 
the  assumptions  of  Lemma  2.  Specifically,  we  investigate  the 
conditions  under  which 


0 


U^A  +  Au  +  /{A}  +  0,  -Ao-A 
mm  12 


A 

R2  >  0  ,  0i  >  0 


y  (86) 


possesses  a  unique  positive  definite  solution  for  which 
Aq  [  (Uj^-02A)  Q+Q  ( Um“C!2  A)  {Q  }]  is  asymptotically  stable. 

Before  proceeding,  however,  it  is  necessary  to  state  various 

preliminary  results.  The  first  of  these  concerns  an  elementary 

(45) 

property  of  positive  semi-def mite  matrices 
Lemma  3 

Every  sequence,  (x^},  of  hermitian  positive  semi -definite 
matrices  bounded  below  (above,  resp. )  with  negative 

semi-definite  (positive  semi-definite,  resp.)  for  each  i  converges 
to  a  positive  semi-definite  limit. 
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T 


The  second  result  concerns  the  preservation  of  detectability 
(Theorem  3.6{ii)  of  Reference  (44)): 


Lemma  4 

k 

If  M  5  0  and  (A,M  )  is  detectable  (reconstructible)  then  for 
all  Q  >  0,  N  >  0  and  any  B  and  F  the  pair  (A  +  BF,{M  +  Q  +  F^NF)^) 
is  detectable  (reconstructible) . 

Further,  it  is  well  to  recall  a  basic  result  for  quadratic 
optimization  of  deterministically  parametered  systems  (Theorem 
12.2  of  Reference  (44)): 


Lemma  5 


Under  previous  definitions  of  and  6#  consider 


A  A  A 

0  =  +  hu  +  S  -Ao-A 

mm  2 


^  A  ^  -l^H 


(87) 


with 


S  2  0  ,  R2  >  0 

A  A 

(IT  ,  B)  stabilizable 
m 

—  ^ 

(u^,  S  )  completely  reconstructible 

Then  (87)  has  a  unique  positive  definite  solution,  A^,  and 
(Uj^-o^A^)  is  asymptotically  stable. 

The  final  preliminary  result  illustrates  the  importance  of 
stabilizability  and  is  indispensable  to  the  main  Theorem  of  this 
section : 
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Lemma  6 

With  3,  I  and  as  previously  defined  and  (p,  B)  stabi- 

lizable,  there  exists  a  Y  >  0  such  that 


A^[A(p  -ijj  -a2Y)  +  (u  -*5i  -a^Y)^h  +  I{A}] 


is  asymptotically  stable.  If,  in  addition,  the  system  possesses 
no  rigid  body  modes,  Y  may  be  chosen  as  any  positive  diagonal 
matrix. 

The  proof  is  given  in  Appendix  6.  With  this  lemma  we  are  in 
a  position  to  demonstrate  the  main  result. 

Theorem  9 

_ 

Assuming  (u  detectable,  adopt  decomposition  (79)  and 

m  1  ^  ^ 

define  s  0,  o  y  and  as  in  Lemma  1  so  that,  in  particular, 
(tm/Oi^)  is  completely  reconstructible.  Then,  if  (p,  g)  is 
stabilizable : 

^  A  A  A 

0  =  +  r{A}  +  -AO2A  (88) 


has  one  and  only  one  positive  definite  solution,  A^»  and: 


A 


-02.)^)  Q  +  Q(u^  +  HQ}] 


and  asymptotically  stable. 


64 


A 


Proof 

To  show  existence  of  a  positive  definite  solution  we  define 
an  infinite  sequence  {A^^}  by: 


A  A  A  A  A 

0  =  T.  +  O-  +  A.U  +  U^A.  +  J{A.}  -A.O-A. 
1  1  1  m  mi  1  1  2  1 


a .  = 
1 


A 

u. 


m 


A 

-a_A 


2  i 


l.a.  +  a^E.  +  HI.) 

11  11  1 


T. 

1 


a. 


(89) 


A 


i+1 


Z  . 

1 


d. 


J 


for  each  i.  First,  assume  there  exists  a  positive  definite  A^ 
for  which  is  positive  semi-definite.  Then,  rearranging  (89a), 
we  have  that  exists  as  a  positive  definite  solution  of 


^A. 

m  1 


A.U  +  s  -A.o_A. 
1  m  1  2  1 


a. 


where 

A  A  A 

s=T.  +a^+7{A^}>0 
—  h 

(u  ,  s^)  completely  reconstructible 

_  A 

(Ujj^,  6)  stabilizable 


b. 


(90) 


c . 

d. 


Note  that  ^  so  that  the  unstable  subspace  of  u 

contains  that  of  IT  and  (90d)  is  implied  by  stabilizability  of 
/V  in 

(IT,  B)  .  Moreover  (90c)  is  implied  by  Lemma  4  (setting 
Q  =  T.  +  r{A.}).  Conditions  (90b, c,d)  conform  to  the  assumptions 

Lx  ^ 

of  Lemma  5;  whence,  asymptotically  stable. 
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A 


Rearranging  (89a)  once  again,  is,  by  assumption,  a  posi¬ 
tive  definite  solution  of 


A  A 

0  =  A^(a  -Jj/)  +  (a  +  /{A^}  +  s 

where ; 

A  ''  ^ 

a  =  u  ~^2^i 

A  A 

S  =  T.  +  +  Aia2A.  2  0 

_  *  w 

(a  -HT I  s  )  completely  reconstructible 


>  (91) 


-  ''  h 

The  last  condition  follows  from  reconstructibility  of  (Um/Cf-,  ) 

—  1 

and  application  of  Lemma  4  (setting  N=R-  >0,F=6A. r 

A.  ^  ^  X  A 

Since,  in  addition,  (a  -HD  =  ^jn"^2^i 


Q  =  T^  5  0  and  B  =  -0^2  ^ 
is  stable.  Theorem  7  implies  that 


^A^^'^i  °‘i^  ^  i(A)] 


is  asymptotically  stable. 

Elimination  of  T^  from  (89a)  and  (89c)  gives: 


0  =  + 


(92) 


H  ^ 

Since  Ay^[Aa^  +  A  +  J{A}]  is  asymptotically  stable,  and 
A^^j^  are  seen  to  be  negative  and  positive  semi-definite, 
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respectively,  by  Theorem  6.  Moreover  reconstructibility  of 

(U  ,0,^)  implies  (a. , (o,  +  A.a„A.  +  reconstructible 

m  i  1  1  1  1  1+1 

so  that: 


^itl  = 


a  ■  t  A 
dt  e  [o 


,  +  A . o„A .  + 
1  1  2  1 


"^"i+1 


ait 

)]e^ 


is  nonsingular.  Manipulation  of  (89)  yields: 


T.  , .  =  E . a_E . 
1+1  1  2  1 


(9 


and  is  therefore  positive  semi-definite. 

Thus,  in  summary,  under  the  assumption  that  A^  >  0  exists 
for  which  T^^  in  (89a)  is  positive  semi-definite,  it  follows  that 
(A^^j^  “^i^  “  ^i+1  ^  "^i+l  "  other  words: 


A  =  lim  A. 
®  .  1 
X-^00 


exists  by  Lemma  3.  Since  T^4-0  this  limit  is  a  positive  definite 
solution  of  (88)  for  which  Aq  ( ( yjjj-02  A^)  Q+Q  03  A^)  7  {Q }  ]  is 

asymptotically  stable. 

It  remains  to  show  that  A  >  0  exists  for  which  T  > 

o  o 

Lemma  6,  a  Y  >  0  exists  such  that 


0 .  By 


is  asymptotically  stable,  since  (u,  §)  is  stabilizable.  Now  a 
W  >  0  can  always  be  found  such  that: 


A 

[W 


A 

+  Yu 


m  m 


A  A 

J{Y}  -Y02Y) ] 


>  0 


(94) 


Also,  let  Z  be  the  unique  positive  semi-definite  solution  of 


O-Ziu^  -CjV) 


(u„  -o,V) 
m  z 


H 


Z  + 


A 

J{z} 


+  w 


(95) 


in  accordance  with  Theorem  6.  Letting: 

A  =  Y  +  Z  >  0 
o 

(89a)  becomes 

T  =  [W  -  (a,  +  Yu  +  U^^Y  +  ^{Y}  -Ya_Y)  1  +  Zo-Z  >  0  (96) 

o  1  m  m  z  z 


Thus  Aq  >  0  and  s  0  and  this  completes  the  proof  of  existence. 

To  prove  uniqueness,  let  and  A2  be  any  two  distinct  posi¬ 
tive  definite  solutions  of  (88).  With  reasoning  analogous  to 

that  following  (90)  we  may  conclude  that  both  u  and 

^  ^  m  2  1 

^m  ~^2^2  asymptotically  stable  by  Lemma  5.  Also,  rearranging 
(88)  into  a  form  analogous  to  (91) ,  it  is  seen  that 
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A.[A(p  -a„A,)  +  (u  -o_A  )  A  +  /{A})  and 
/V '  '  m  2  1  m  2  1 

^  ^  ( A I  ]  are  asymptotically  stable 

by  Theorem  7.  Define: 


Z  =  ^2 


Manipulation  of  (88)  yields: 


0  =  Z(m^  -o^A^)  +  (u^  -a,A,)"z  +  J{Z}  -Zo,Z 
m  2  1  m  2  1  2 


0  =  Z(m  -a.A_)  +  (m  -a^h^)  Z  +  I{Z}  tZo^Z 
m  2  2  m  2  2  2 


Since  Za2Z  >  0,  Z  is  both  positive  and  negative  semi-definite  by 
Theorem  6  and  hence  Z  =  0.  □ 

3 . 5  Existence  and  Uniqueness  of  Steady-State  Solutions  - 
Closed-Loop  Stability 

With  the  foregoing  developments  we  are  now  ready  to  determine 
sufficient  conditions  for  existence  and  uniqueness  of  steady  state 
solutions  of  (64) .  The  preceding  results  may  be  combined  and  sum¬ 
marized  as  follows; 

Theorem  10 

Consider : 


2^  _ _  _  _  _ 

-p  =  +  P  Uj^  +  I{p}  +  -pa2P 


p{t^)  =  0 
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where : 


-  A  - 


Oj  i  6R-1b“ 


and  where  M,  1,  6,  and  are  as  defined  in  Section  3.1. 


If: 

"1 

2  0  ,  R2  >  0 

a. 

(U, 

B)  stabilizable 

b. 

(u. 

I5 

)  detectable 

c. 

then: 

A. 


(99)  has  a  unique  positive  semi-definite  solution  for  all 
t  <  t. 


XiXM  — 

B.  ^  P (t)  =  A  where  A  is  the  unique  positive  semi-definite 
solution  of 


0  =  m!!a  +  Au  +  J{A}  +  0,  -Ao„A 
mm  12 


“1  H 

C.  With  K  =  R2  p  r  the  closed  loop  system  is  second  moment 

exponentially  stable  and  almost  surely  exponentially  stable 
for  all  I  >  0. 


(100) 


(101) 


(102) 
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Proof 


Part  A  has  been  shown  in  Lenmia  1,  and  so  we  pass  immediately 

to  Part  B.  First  note  that  ^  so  that  the  unstable 

R  m  R 

subspace  of  is  contained  in  that  of  IT.  Then,  condition  (lOlc) 
implies : 


(u  ,  a^  )  detectable 
m  1 


Thus  we  can  employ  the  reconstructibility  canonical  form  of  the 
system  arising  from  decomposition  (79).  By  Lemma  1,  the  positive 
semi-definite  solution  of  (99)  can  be  reduced  to  the  form 


with  p  satisfying  (82)  ,  where  (y^,  o^^)  is  completely  reconstruct- 
ible.  Since  the  reconstructible  subspace  contains  the  unstable 
subspace  and  (^,  6)  is  stabilizable,  (y,  g)  is  also  stabilizable. 
Thus  the  conditions  of  Theorem  9  are  met  and  (84)  has  a  unique 
positive  definite  solution,  A  ,  and 


is  asympototically  stable.  Then,  by  Lemma  2; 


lim  p(t)  =  A 
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Therefore,  under  conditions  (101),  we  conclude  that: 


lim  p(t) 

t^  too 


(103) 


where  is  the  unique  positive  definite  solution  of  (84). 


Obviously,  (103)  is  also  a  positive  semi-definite  solution 
to  (102) ,  but  it  remains  to  show  that  (103)  is  the  only  such 
solution.  Partition  all  matrices  in  (102)  in  accordance  with 
(79);  in  particular  set: 


A 


A 


(104) 


and  use  (80)  and  (83).  Then  the  (2,2)  sub-block  of  (102)  may 
be  written: 


0  =  u”  A_  +  A_p”  +  I  {A.,} 
mu  2  2  mu  u  2 


-(b“a  +  B^A  )^  R~^(B^A  +  B^A  ) 

12  u  2  2  12  u  2 


y  (105) 


Note  that  if  A  is  positive  semi-definite,  so  is  A 


2‘ 


On  the  other 


«  -  H  ^ 
^2^mu 


hand,  it  was  shown  in  the  proof  of  Theorem  8  that 

'^Q2t^muQ2  Q2^m"u  ^  ^u^Q2^1  '  equivalently,  A^2  f^mu^2  ^ 

[  {Ap}]  is  asymptotically  stable.  Then  from  (105),  Theorem  6  and 

R2  >  0,  A2  is  negative  semi-definite  unless  ^  ^u  2^  van¬ 


ishes,  in  which  case  A^  also  vanishes, 
negative,  we  must  have: 


Thus,  for 


A^  to  be  non- 
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■'2  =  0 


."A,,  =  0 


a. 


:io6) 


Under  these  conditions,  the  (1,2)  sub-block  of  (102)  assumes 
the  form: 


^12  ^ 


At  tU 

12  mu 


whence  Aj^2  vanishes.  Therefore,  any  positive  semi-definite  solu¬ 
tion  of  (102)  must  assume  the  form: 


Lo  0  J 

But  with  this  result,  (102)  reduces  to  (84)  and  (84)  has  been 
shown  to  possess  the  unique  positive  definite  solution  A^.  Con¬ 
sequently  (103)  is  the  unique  positive  semi-definite  solution  of 
(102),  which  completes  the  proof  of  Part  B. 

Because  of  the  result  of  Theorem  9,  the  assumptions  of 
Lemma  2  are  fulfilled-  Thus,  with  the  steady  state  control  gain 
K  =  S^A, 

is  asymptotically  stable  by  Therorem  8  and  the  conclusion  of 
Part  C  follows.  [! 
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The  situation  with  respect  to  the  existence  and  uniqueness  of 
steady  state  solutions  to  (99)  is  thus  seen  to  be  entirely  analo¬ 
gous  to  that  of  the  ordinary  Riccati  equation  for  a  deterministic 
plant.  Note  that  if  the  system  lacks  rigid  body  modes,  stabiliza- 
bility  and  detectability  are  assured  so  that  (101a)  form  the  only 
ancillary  conditions.  In  the  general  case,  we  require  only  that 
the  rigid  body  modes  be  controllable  and  reconstructible  -  condi¬ 
tions  readily  imposed  in  practice. 

In  place  of  guaranteed  stability  for  the  nominal  system  which 
results  from  the  solution  of  the  deterministic  plant  regulation 
problem  we  here  obtain  the  stronger  assurance  of  second  moment 
stability.  This  is  indirectly  a  consequence  of  the  use  of  a  mean- 
quadratic  performance  criterion,  i.e.,  existence  of  a  mean-square 
optimal  control  necessarily  demands  second  moment  stability. 

Furthermore,  due  to  the  particular  form  of  the  uncertainties 
considered,  even  stronger  guarantees  of  stochastic  stability  may 
be  obtained.  The  following  result: 

Theorem  11 

Assume  the  conditions  of  Theorem  10  and  that  ^  is  the  posi¬ 
tive  semi-definite  solution  to  (102).  Define: 


S  =  -^Ip  +  /{pl  -^P/ 


(107) 


Then : 


0  <  S  < 

(u,  S^)  detectable 


(108) 


the  proof  of  which  is  given  in  Appendix  7,  leads  to  the  conclu¬ 
sion  : 

Theorem  12 

With  p  the  positive  definite  solution  to  (102)  under  the 
conditions  of  Theorem  10,  then  (IT  -a^P^)  is  asymptotically  stable 
and  the  closed  loop  stochastic  system  is  almost  surely  exponen¬ 
tially  stable  and  p*"^  mean  exponentially  stable  for  all  p  >  0 
and  all  values  of  the  decorrelation  times. 

Proof 

By  rearrangement  of  (102) ,  p  is  seen  to  exist  as  the  posi¬ 
tive  semi-definite  solution  of 


0  = 


-H- 

U  P 


+  p  M  + 


S-0O2P 


(109) 


_  _  p 

where  (u ,  6)  is  stabilizable  by  assumption  and  (u»  S ^)  is 
detectable  by  Theorem  11.  Consequently  (see  Theorem  12.2  of 
Reference  (44))  (y  asymptotically  stable.  Since,  in 

addition,  the  random  portion  of  y  is  sJcew-hermitian  (see 
Appendix  6)  the  remaining  conclusions  follow  by  application  of 
Criterion  6  of  Reference  (46) .  □ 

Thus,  the  stochastic  regulator  design  guarantees  a  very 
considerable  measure  of  stochastic  stability. 

3 . 6  Constant  Gain  Stochastic  Design;  Computational 
Procedures 

Because  of  the  relative  ease  with  which  constant  gain  con¬ 
trollers  may  be  implemented,  attention  is  henceforth  restricted 
to  the  steady  state  case.  In  the  following,  p  denotes  the 
positive  semi-definite  solution  to  the  steady  state  version  of 
(99)  : 
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0 


-H- 


I{p}  +  0^ 


—  Li  —  1  T 

=  U  -HI 


02  =  BR2  B 


(110) 


under  the  conditions  of  Theorem  10.  In  this  section  various  con¬ 
vergent  iterative  procedures  for  the  numerical  solution  of  (110) 
are  set  forth. 

The  first  method  is  suggested  by  Theorems  9  and  10: 

Theorem  13 

Using  the  decomposition  of  Lemma  1,  choose  Y  >  0  such  that: 


A  A  A  A 


AA[A(Mm  -^2^^  ^ 


is  asymptotically  stable  and  define  W  >  0  and  Z  >  0  so  that: 


A  A 


-H 


[W  -(O^  +  +  I{Y}  -Ya2Y)  ]  >  0 


(111) 


A  A  A  A 

0  =  Z(i7  -o^.)  +  (m  -o_Y)”z  +  I{z}  +  w 
m  ^  m  z 


(112) 


Then : 


lim  =  p  (113) 

kf  00 
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where : 


“  *  Vl  ‘>k+l“'n,  ■°2''k> 


*  "‘“k+l*  ■"  °1  *  V2Pk 


P 

o 


Y+Z  0 
0  0 


y  (114) 


Proof 

With  (80)  through  (83)  from  Lemma  1,  induction  on  k  from  (114) 
shows  that 


i- 


for  all  k  >  0.  The  sequence  (Pj^}  is  seen  to  be  defined  by  rela¬ 
tions  (89) .  Furthermore,  (111)  ,  (112)  and  the  choice  for  Y 
ensure  that  defined  by  (96)  is  positive  semi-definite.  Refer¬ 
ring  to  the  proof  of  Theorem  9,  this  guarantees  the  convergence 
of  the  sequence  (89).  Thus,  sequence  (114)  converges  to  a  posi¬ 
tive  semi-definite  limit.  By  (103)  this  limit  is  the  solution 
of  (110).  □ 

The  above  result  corresponds  to  the  Newton-Raphson  method  and 
therefore  exhibits  quadratic  convergence  in  a  sufficiently  small 
neighborhood  of  the  solution  of  (110) .  The  primary  difficulty 
arises  from  solution  of  a  stochastic  Lyapunov  equation  on  each 
iteration.  Also,  of  course,  proper  choice  of  the  starting  value 
may  occasion  some  inconvenience. 
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i 

Ambiguity  in  the  starting  value  is  removed  in  the  following 
method. 

Theorem  14 


lim 
kf  “ 


Pk 


P 


(115) 


where  the  positive  semi-definite  sequence 


(Pk) 


is  defined  by 


°  '^m^k+l  ^k+l^m  °1  ■'^k+l°2Pk+l 


Po  =  0 


(116) 


Proof 


Using  the  representation  of  Lemma  1,  (116)  is  seen  to  yield: 


Pk  = 


\  0 


(117) 


for  all  k  >  0.  By  substitution,  the  are  found  to  satisfy: 


A  ^  A  A  A 

0  =  m"A,  ,1  +  A,  .iM  +  J{A,  }  +  o,  -A,  ..OtA,  , 

m  k+1  k+l^m  k  1  k+1  2  k+1 


A  =  0;A,  >0,Vk>0 

o  k 


(118) 
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Now  let: 


P  = 


0 


0 

0 


be  the  positive  semi-definite  solution  to  (110)  where  A  is  the 
positive  definite  solution  of  (88)  and  asymptotically 

stable  by  Theorem  9.  (88)  and  (116)  may  be  manipulated  to  yield: 


0  =  (u^  -o^A)  t  (A-A^,i)(u^  -o^A) 


+  r{A-Aj^}  +  (  A-Aj^_^^)  02  ( 


(119' 


Thus,  by  Lemma  12.1  of  Reference  (44),  if  A-A^^  is  positive  semi- 
definite  then  so  is  A-Aj^_^^.  Furthermore  if  Aj^  is  positive  defi- 

'1 

Thus  Lemma  5  is  applicable  to 


nite,  then  because  (U^,  o^'^)  is  completely  reconstructible ,  sc  is 


m 


(Pm'  (^i  +  Lemma  4. 


(118)  ,  whence  ”^2^)c+l  asymptotically  stable.  Manipulation 
of  (118)  yields: 


^  ^^m  "°2^k+l^  ^^k+l~^k^  ^^k+l"^k^  “°2^k+l^ 


^^^k  ^k-1^  ^^k+l"^k^®2^^k+l  ^k^ 


(120) 
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so  that  if  (Aj^  positive  semi-definite  then 

is  also  positive  semi-definite.  Thus,  by  Lemma  3,  the  sequence 
’Aj^}  converges  to  A  provided  that  -A^  >  0.  But  this  is  assured 
by  the  choice  =  0  ^nd  by  virtue  of  Lemma  5.  Consequently,  with 
(117),  the  sequence  converges  to  the  positive  semi-definite 

solution  of  (110) .  □ 

This  method  may  be  termed  the  "perturbed  weighting"  method 
since,  in  effect,  it  involves  sequential  modification  of  the  state 
weighting  matrix,  o^.  Obviously,  the  problem  of  determining  a 
suitable  starting  value  is  eliminated.  Moreover,  each  iteration 
involves  solution  of  the  standard  Riccati  equation.  On  the  other 
hand,  this  method  exhibits  only  linear  convergence  in  contrast  to 
the  quadratic  convergence  of  the  Newton  method. 

3 . 7  Asymptotic  Properties  for  Large  Uncertainties 

Having  determined  that  a  unique  positive  semi-definite  solu¬ 
tion  of  (110)  exists  for  all  positive  I,  it  is  natural  to  enquire 
what  behavior  "p  attains  for  large  uncertainties,  i.e.,  for  very 
small  decorrelation  times.  In  examining  this  question  we  are 
mainly  concerned  with  the  resulting  form  of  the  control  law  for 
the  elastic  modes.  Thus,  in  the  remainder  of  this  report,  it  is 
assumed  that  the  system  possesses  no  rigid  body  modes,  i.e.,  that 


u  =  ih 


n  =  n^.  This  restriction  permits  a  relatively  straight¬ 


forward  development  and,  moreover,  the  qualitative  form  of  the 
results  to  be  derived  below  is  expected  to  hold  even  in  the 
presence  of  rigid  body  modes. 

The  principal  effects  of  frequency  uncertainty  for  small 
values  of  the  decorrelation  times  are  connected  with  the  distinct 
character  (noted  after  the  proof  of  Theorem  4)  of  the  diagonal 
and  off-diagonal  portions  of  equation  (110) .  Since  we  mean  to 
treat  the  diagonal  and  off-diagonal  elements  of  p  separately,  it 
is  expeditious  to  introduce  the  notation: 
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<M>  £  M  -{M}  (121) 


for  any  square  matrix  M.  In  addition,  define: 


{p*}  may  be  recognized  as  the  solution  for  the  diagonal  portion 
of  ^  obtained  from  the  diagonal  portion  of  (110)  with  ignored. 

With  these  definitions,  we  first  establish  certain  bounds  on 
and  as  follows; 

Theorem  15 

Under  the  conditions  of  Theorem  10,  and  with  "p  the  positive 

semi-definite  solution  of  (110)  with  11  =  11: 

e 


l<P>kil®  _i  •_ - ^ - 


B.  Defining  Bj^  _  ^  Kp>02} 


(123) 


=  Kp>0  2<P>)]<. 


(124) 


'"I’k  -Cfc  "  0 


(125) 
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then : 


I  B^+2{p*},,{a2}]^-2Reu 


C. 


If  (o2^kk  ^  then: 


i(p}fc-(p*),^ 


<  2 


{\ 


Ikk 

^2kk 


[r.  +  r; 


«,j^k  I 


Proof 

As  a  preliminary  step,  decompose  each  term  of  (110)  into 
diagonal  and  off-diagonal  elements  to  obtain: 


0  =  2Reyj^{p}j^  + 

0  =  <P>kjtM*+yj  -Jsdj^+ij)]  +  <Pi>kj-<PP2P>kj 


Now  consider  part  A.  Since  c ^  >  0 ,  p  o 2P  is  positive  semi 
definite.  Consequently  {F  non-negative  and: 


l<pa2P>kjl 


—  —  _  —  k 

<  ({pa2P}j^{pa2P}  j) 


k  5^  j 


(126) 


(127) 


(128) 


its 


(129) 


(130a) 


Also,  from  (129a): 


{po2P}j^ 


since  Reu  <  0.  Then  {130a)  yields: 


KP02P\j 


Using  this  inequality  in  (129b)  and  noting  that  because  > 


'Ikj 


we  obtain  (123) . 

With  the  above  results,  inequalities  on  may  be  shown, 
starting  with  part  B.  Expand  the  quadratic  term  in  (129a) 
to  get: 


with  and  Cj^  as  defined  in  (124).  Defining: 


A  =  {p}  -{p*} 


(130b) 


(131) 


(132) 


(133) 
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and  using  (132)  with  (122),  we  obtain  for  A: 


[Bj^+2{p*}^{02},^-2Reyj^]Aj^  +  Cj^+{p*>kBj^  =  0 


(134) 


Solving  (132)  for  and  (134)  for  Aj^: 

(p), . 

J?  [-a±/a^-4B]  ^  0 

a. 

\  = 

!3[-a±/a^-48] 

b. 

a 

a. 

a 

=  a  +  2{p*}j^ 

b. 

8 

°  ‘V‘'=*'k®kl 

c. 

(135) 


(136) 


where  the  same  sign  appears  before  the  radical  in  both  (135a) 
and  (135b) .  Suppose  that  the  negative  sign  is  to  be  taken  in 
(135).  Then,  (135a)  implies: 

2{p*}j^  >  a  +  /a^-A  6 


or;  rearranging: 


0  >  a  +  /ot^- 


4[V{ai}J/{a2}k 
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But  in  view  of  condition  (125) ,  the  right  hand  side  above  is 
inherently  positive,  implying  a  contradiction.  Therefore,  the 
positive  signs  must  be  chosen  in  (135) . 


a 

B 


[-a  + 


{a2^k 


^°2^k 


/a^-AB  ] 

[Bj^-2ReUj^]  + 


2{p*}^ 
>  0 


a. 


b. 


(137) 


c. 


Note  that  (132)  may  be  rewritten: 


0  =  (^-{02}  +  {p}  (y-{a2}{p}-{a2<P>}) 

+  {■p}^{a2}  +  {a^}-{<F>a2<^>} 


Now,  {p}  exists  as  a  positive-definite  matrix,  while 
{o^}  -{<p>a2<P>}  is  positive  by  assumption  (125) .  Consequently 
(u  -{02H^}  asymptotically  stable,  i.e.,  a  is  pos¬ 

itive. 

With  this  property,  the  inequalities: 

1  -  <  1  +  ^  ,  V  X  real 
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may  be  used  in  conjunction  with  (137)  to  obtain: 


A 


k 


iBi 


(138) 


With  substitution  of  (137b, c),  the  result  (126)  follows. 
Part  C.  can  be  shown  more  easily.  From  (132) : 


{a2}j^{p}j^  2  l^k^M  '*■  2ReMj^] 

{a2}ktp}k  -  ^^i^k  2Repj^] 


^  (139) 


where : 


^^k^M  =  ^  «-7k  ’^^2kk'^2J!,«,' P«,k'M 

’'"^Tu^imml  ^ilkUl  ^mkU 


(140) 


and  denotes  the  upper  bound  given  by  (123).  Note  also 

that : 


Consequently,  the  second  inequality  of  (139)  actually  gives  an 
upper  bound  for  From  {139b): 


\lM'*‘2R»l'k)^+4(Ol)k<»2>k> 


SO  that: 


{P},-{PM,^| 


2{^ 


[IB 


k'M 


Vl 


B, 


|j^+4{0i)j^{02}j. 


-  ^ 


°l>k^"2>lcJ  (141) 


Using : 


/I+x  -1  <  I  ,  ¥  X  >  0 


in  (141)  and  substituting  expressions  (140)  ,  the  result  (127)  is 
obtained.  □ 

Of  the  bounds  given  in  Parts  B.  and  c.  ,  C.  is  the  simpler 


but  does  not  apply  for  cr2j^j^  =  0. 


On  the  other  hand  (126)  is 


applicable  for  sufficiently  high  uncertainty  levels  and  a 


2kk 


>  0 


and  gives  a  closer  bound  for  small  values  of  a 


2kk' 


The  presence  of  (I,  +1.)  in  the  denominator  of  (123)  clearly 

^  3 

shows  that  one  effect  of  frequency  uncertainties  is  Iro  suppress 
the  off-diagonal  elements  of  the  expected  cost.  In  this  connec¬ 
tion,  use  of  the  above  bounds  easily  leads  to  the  following  con¬ 
clusions  : 
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Theorem  16 


88 


It  is  clear  that  the  bound  given  by  (123)  decreases  mono- 
tonically  with  increasing  J.  Hence: 


lim  |<p  ^  >1  =  0  (146) 


Moreover,  there  exists  a  J  sufficiently  large  that  condi¬ 
tion  (125)  is  satisfied  and  the  bound  (126)  may  be  em¬ 
ployed.  In  view  of  (146)  ,  the  limit  of  the  right  side  of 
(126)  is  zero.  □ 

From  the  definition  (lOe)  of  and  the  form  of  <!>  given  by 
(8) ,  it  is  seen  that: 

^'^l^  2m-l ,  2m-l  ~^*^1^2m,2m;m=l,...,n  (147a) 


Similarly,  from  (10b),  (8)  and  (49): 
^®2^  2m-l,2n-l  ^  ^'^2^2m,2m  ' 


(147b) 


Consequently,  from  (12b)  and  (122): 

^P*^2m-1  "  ^P*^2m  ' 


and  K  may  be  written: 


X  A0-A099  III  MASSACHUSETTS  ZNST  OF  TECH  LEXINSTOH  LINCOLN  LAB  F/6  13/lS 

OPTIMAL  RE6ULAT10N  OF  STRUCTURAL  SYSTEMS  MlTH  UNCERTAIN  PARAMET-^CfUl 
FEB  ai  DC  HYLAND  Fi9628-80-C-0002 

UNCLASSIFIED  TR-551  ESD-TR-80-246  NL 


ESD-TR-80-246 


K 


—  It 

R-  B  X  block-diag 


2oj2  0 
in 


0 

1 


(148) 


Thus,  since  the  odd  rows  of  B  vanish,  the  odd  columns  of 
K  also  vanish.  □ 

e.  Define: 


=  >5c“{p*}?  >  0  ,  V  ^  ji  0 


Then,  for  the  system  (145): 

=  ({p*}^5)“[Rey  -{p*}^BR2^B”{p*}^]{p*}’*C 
<  0 


since  y  <  0  and  R2,  {p*}  >  0.  Therefore,  JI  is  a  Lyapunov 
function  for  (145)  and  the  stated  conclusion  follows.  D 

The  stochastic  Riccati  equation  does  indeed  have  an  asymptotic 
solution,  of  a  particularly  simple  form,  in  the  limit  as  all  decor¬ 
relation  times  tend  to  zero.  According  to  Parts  B  and  C  of  the 
above  Theorem,  the  resulting  rate  feedback  control  law  renders  the 
closed-loop  system  stable  for  all  values  of  the  open  loop  fre¬ 
quencies.  This  feature  appears  to  be  a  natural  consequence  of  the 
inclusion  of  pareuneter  uncertainties  as  an  intrinsic  part  of  the 
model  and  illustrates  the  qualitative  principle:  "If  nothing  is 
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known  regarding  the  open  loop  modal  frequencies  then  the  mean- 
square  optimal  control  is  the  control  law  which  is  inherently 
dissipative  under  frequency  uncertainties."  Furthermore,  Part  C. 
shows  that  the  asymptotic  form  of  the  control  also  guarantees 
stability  in  the  face  of  uncertainties  in  the  modal  damping  ratios. 
Thus  a  greater  degree  of  robustness  has  been  obtained  than  was 
originally  sought. 

3. 8  Asymptotic  Properties  for  High-Order  Modes ; 

Incoherence  and  Decoupling 

Theorem  16  pertains  only  to  the  case  in  which  the  decorrela¬ 
tion  times  are  all  uniformly  small.  Analogous  results  are  to  be 
expected  when  uncertainties  in  the  low  frequency  modes  are  small 
while  modelling  accuracy  degenerates  for  modes  of  increasing 
order.  Supposing  the  modes  to  be  arranged  in  order  of  increas¬ 
ing  nominal  frequency,  we  anticipate  that  correlation  between 
distinct  modes  will  be  suppressed  to  a  greater  and  greater 
degree  the  higher  are  the  orders  of  the  modes  involved.  Simi¬ 
larly,  the  expected  cost  matrix  will  become  increasingly  diago¬ 
nalized  toward  the  lower  right  hand  corner  and  that  portion  cor¬ 
responding  to  very  high  order  and  poorly  known  modes  will  approach 
the  asymptotic  form  (143) . 

To  assess  the  degree  to  which  off-diagonal  elements  of  P  are 
reduced  in  magnitude,  these  must  be  scaled  properly  relative  to 
the  corresponding  diagonal  elements.  Thus,  for  the  purpose  of 
the  following  development,  we  introduce  the  co-state  correlation 
coefficient  matrix,  p: 

p  t  {F}  **  ^{p}  ^ 
so  thati 

“  ^kj^’^^kk^j  j 


a. 


b. 


y  (149) 
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is  clearly  analogous  to  the  state  correlation  coefficient  formed 
from  the  covariance  and  has  similar  properties. 


With  (149)  ,  and  as  a  consequence  of  the  bounds  given  in 
Theorem  15,  one  obtains: 


Theorem  17 

Arrange  modes  in  order  of  increasing  nominal  frequency  and 
suppose  that  n,  the  number  of  modes  retained  in  the  model,  is 
arbitrarily  large.  Considering  only  the  reconstructible  modes 
partition  the  state  vector  thus: 


C  = 


2N, 


(150) 


and  define  the  corresponding  partitions  of  p  and  ^  by: 


■N 


(151) 


j 


1 


Assume  that  the  damping  coefficients,  rij^  are  non-zero  and 
bounded  for  all  )c;  that  the  input  matrix,  B,  and  the  control 
weighting  matrix,  R2,  are  bounded  from  above  and  from  below, 
respectively,  and  the  following  restrictions: 


92 


1.  monotone  increasing  with  increasing  k 
-1  ‘^N(k)  ^  '^Ikk  "  ‘^N(k)  '  ^ 


where  >  0 


1  «- 


3.  —  ““m)  nondecreasing  with  m 


Then,  given  e  >  0,  there  exists  an  sufficiently  large  that: 


i<Pi>jki  <  p 


"^'jk' 


<  e 


b. 


c. 


>  (152) 


where  {p|}  is  (143)  evaluated  for  the  "I"  modes.  Moreover,  each 
of  the  quantities  on  the  left  of  (152)  are  bounded  by  monotoni- 
cally  decreasing  functions  of  k . 


The  proof  is  given  in  Appendix  8. 


Conditions  (1)  through  (3)  of  the  above  Theorem  demand  par¬ 
ticular  note.  Clearly,  condition  (1)  demands  that  the  decorrela¬ 
tion  times  expressed  as  multiples  of  the  corresponding  natural 
periods  of  vibration  decrease  with  increasing  modal  order  - 
reflecting  a  general  decline  of  modelling  accuracy  for  the  higher 
order  modes.  Condition  (2)  assumes  an  "energy  weighting"  on  the 
state  and  represents  a  rather  extreme  assignment  of  cost  for  the 
high  order  modes.  Often  in  practice,  displacements  alone  would 
be  included  in  the  performance  index,  resulting  in  a  more  modest 


rate  of  increase  of  with  k.  Finally,  condition  (3)  requires 

that  the  modal  frequency  separation  considered  as  a  function  of 
frequency  be  bounded  from  below  by  a  constant,  i.e.,  that  modal 
density  decrease  with  increasing  frequency.  This  condition  is 
also  likely  to  hold  in  practice,  at  least  above  a  certain  frequency. 


Inspection  of  Appendix  8  reveals  that  the  rates  of  increase  of 
modal  density  and  modal  state  weighting  for  which  Theorem  17 
retains  validity  are  related.  In  particular,  for  increasing 

as  a  power  of  less  than  two,  the  above  results  will  hold  and 
will  retain  validity  even  when  modal  density  increases  with 
increasing  frequency.  In  any  case,  however,  the  Theorem  as 
stated  is  applicable  to  most  situations  of  interest. 


We  may  say,  in  view  of  (152) ,  that  the  partitioning  (150) 
apportions  the  modes  into  "coherent"  and  "incoherent"  systems. 
This  assignment  of  modes  is  uniquely  defined  by  the  maximum 
tolerable  correlation  level,  e,  which  may  typically  be  set  <<1. 
Further,  may  be  termed  the  "coherence  limit"  associated  with 
the  correlation  level  e.  Thus  the  qualitative  features  of  the 
mean-square  optimally  controlled  system  are  analogous  to  those 
identified  in  Section  2.2  for  the  uncontrolled  system. 


The  coherence  limit,  N^{e),  (where  we  now  display  the  func¬ 
tional  dependence  upon  e  explicitly)  may  be  estimated  through 
use  of  the  bounds  given  by  Theorem  15.  However,  useful  and  con¬ 
venient  results  may  be  obtained  with  approximate  forms  of  these 

bounds.  Very  simply,  if  e  is  small,  then  for  k  ^  2N  ,  {^},  may 

V,  .  c  K 

This 


be  approximated  by  in  evaluating  Kp\j  I  from  (123). 

gives  the  following  approximate  determination  of  N^(e): 
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_ I 


N^(e);  |<p>. 


2Nc.j 


<  _ c _ _ _ c _ -* 

^  -‘‘I^2N  +  ^j)I 

C  c 


<  £  ,  V  j  (153) 


Furthermore  if  for  k  ~  2N^  the  damping  is  known  to  be  small  in 
the  sense: 


(154) 


we  have  {p*}]^  —  (^ikk'^°2kk^  that  (153)  becomes 


N^(e\:  !<P>2N^,j 


2  (f<^1^2N^^‘^2^2N^^‘^l^  j^°2^  * 


<  £  ,  V  -j  (155) 


The  above  expressions  (particularly  (155))  reveal  the  depend¬ 
ence  of  N^(e)  upon  the  various  design  parcuneters.  Clearly  if 
higher  control  authority  is  demanded  (i.e.,  {02^k  ^'^l^k 

creased)  N„(e)  is  increased.  Most  important,  N  (e)  strongly 
c  c 

depends  on  the  rate  of  increase  of  the  reciprocal  decorrelation 
times,  Ij^,  with  k,  and  the  more  rapidly  increasing  are  the 
the  smaller  is  N^(e).  In  fact,  the  dependence  of  1^,  upon  k  might 
be  chosen  in  an  ad  hoc  fashion  so  as  to  achieve  a  desired  coher¬ 
ence  limit. 


95 


Obviously  if  e  is  sufficiently  small,  is  negligible  so 

that  modal  frequency  uncertainties  effectively  decouple  the 
coherent  and  incoherent  systems.  Moreover,  p'^  is  approximately 
{p*}  as  given  by  (143).  Thus,  for  the  incoherent  system,  (110) 
naturally  gives  rise  to  a  rate  feedback  control  law  which  is 
stable  regardless  of  the  values  of  modal  frequencies  or  damping 
ratios.  For  sufficiently  large  values  of  the  reciprocal  decor¬ 
relation  times  and  fixed  e<<l,  N^(e)  dirinishes  to  less  than 
unity;  the  whole  system  becomes  incoherent  and  the  entire  control 
law  reduces  to  the  asymptotic  form  given  in  Theorem  16. 

The  above  properties  can  greatly  simplify  computation  of  the 
control  gain  for  system  models  of  very  high  dimension.  We  first 
note,  without  dwelling  on  the  matter  in  any  detail  here,  that  the 
expressions  given  in  Theorem  15  or  their  approximate  counterparts 
may  be  used  to  bound  the  contributions  of  and  “{pj}  to  k. 
These  contributions  decline  with  increasing  and  are  convergent 
even  as  n  tends  to  infinity. 

Thus,  under  the  assumptions  of  Theorem  17,  the  calculation  of 

<  may  proceed  as  follows.  First,  using  the  bounds  given  in 

Theorem  15,  determine  N  sufficiently  large  that  the  magnitude  of 

c 

all  elements  in  the  error; 


(156) 


incurred  in  the  approximation; 


{pp 


(157) 


may  be  considered  negligible.  Then  compute  <  from  (157)  ,  deter¬ 
mining  P'^  from 


0  =  D  + 
cm  c 


p  y  +  I  {p  }  + 
c  cm  c  c 


o,  -p  a_  p 
Ic  ^c  2c^c 


(158) 


by  either  of  the  methods  of  Section  3.6  and  calculating  {p*}  from 
(122) . 

(158)  is  the  upper  left  sub-block  of  (110)  when  partitioned 

in  accordance  with  (150) ,  and  is  a  stochastic  Riccati  equation  of 

order  2N„  2N  .  If  the  number  of  coherent  modes,  N^,  is  modest, 
c  c  c 

either  the  Newton  method  or  the  perturbed  weighting  method  may  be 
employed  at  reasonable  computational  cost.  The  calculation  of 
{p£},  since  it  proceeds  from  simple  analytical  expressions,  entails 
no  computational  burden  whatsoever.  Note  that  the  error  incurred 
in  Pj  —  {pj}  is  greatest  on  the  upper  left  corner  and  diminishes 
toward  the  lower  right.  Consequently,  somewhat  above  the  coherence 
limit,  as  many  modes  as  desired  may  be  accommodated  without  incur¬ 
ring  significant  computational  effort. 

We  conclude  that  if  the  modeled  uncertainties  increase  with 
modal  order  and  the  number  of  coherent,  relatively  well  known  modes 
is  modest,  the  solution  of  the  stochastic  Riccati  equation  and  the 
control  gain  may  be  computed  with  satisfactory  accuracy  for  systems 
of  arbitrary  order. 
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3. 9  Efficient  Computational  Procedures  for  the  Coherent 
System 

Let  us  retain  the  conditions  of  Theorem  17  and  consider  the 
stochastic  Riccati  equation  for  the  coherent  system: 

®  ~  ^cm^c  ^  ^c^cm  ^  ^  ^Ic  ^c®2c^c 

P(,  e  C 

Here  we  explore  the  numerical  treatment  of  (158)  when  the  number 
of  coherent  modes  is  too  large  to  permit  efficient  application  of 
the  methods  of  Section  3.6. 

In  this  connection,  the  diagonalization  effect  of  frequency 
uncertainties  noted  previously  suggests  a  novel  iterative  numeri¬ 
cal  technique; 


—  d  —  —  _ 


Po  =  {P*} 


r  (159) 


where,  for  convenience,  the  "c"  subscript  has  been  dropped.  The 
rationale  for  (159)  can  be  understood  through  examination  of  the 
diagonal  and  off-diagonal  portions  separately; 
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0 


0 


(U  -HI  -to,Pv^) 


k+r 


2'^k 


*  <°1>  +  ‘'’2'‘k*“<''k>  <“k>‘'’2''k>  -<V2‘’k> 

<“  -‘»2“k>>"‘^k+l'  "■  ‘^k-n’*^  -*°2Pk’> 


+  {Oj^}  +  2Re{pj^}  {a2Pj^}  ■{Pk02'^k^ 


(160) 


where  ^P©^  ~  {p*} 


<»<,>■  » 


y  (161) 


At  least  for  modes  "corresponding  to  the  lower  right  entries  of  'p , 
the  reciprocal  decoi  ilation  times  will  be  relatively  large.  Thus 
the  presence  of  -HI  --u  the  coefficient  matrix  of  in  (160a) 

suggests  that  <  will  be  "small"  in  the  sense 


l<Pk+i>l 


<  < 


Note  that  (y  -HI  “^^P2^k^^  eigen¬ 

value  matrix  of  the  mean  system  evaluated  to  within  a  first  order 
perturbation  in  the  control.  Equation  (160a)  can  thus  be  seen  as 
a  perturbation  approximation  to<^>for  the  mean  system  with  its 
decorrelation  damping  [-HI) . 

V  _  _ 

If  indeed  j  |<P>1  1  is  small,  the  contribution  of  ^  to  Pj^  in 
(160b)  will  also  be  "small"  so  that  provides  a  perturbation 

on  what  would  be  Newton's  method  for  determination  of  in  the 

case  <p>=  0.  Obviously,  to  within  the  approximation  that  <  pj^>  is 
negligible,  (160)  yields  =  {^  }  =  {"p}  since  {p*}  satisfies 

(158)  identically  in  the  case<^>=  0. 


The  above  reasoning  motivates  our  use  of  the  term  "asymptotic 
refinement"  for  the  method  of  (159)  since  it  represents  a  sequen¬ 
tial  correction,  by  simple  substitutions,  of  the  asymptotic  solu¬ 
tion  given  by  Theorem  16. 
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One  clear  advantage  of  (159)  is  that  the  coefficient  matrices 
of  ^  Lyapunov  equations  which  must  be 

solved  on  each  iteration  are  already  diagonal.  Thus,  each  ele¬ 
ment  of  is  determined  separately: 


Mj 


(162) 


{Ol}^t2Re{p^^^Ja^p^}^-{p^o^p^^} 

-2Re(yj^+{02^j^}) 


Thus,  when  it  does  converge,  the  asymptotic  refinement  method 
is  considerably  faster  than  the  methods  of  Section  3.6  since, 
apart  from  evaluations  of  the  form  (162)  it  requires  only  a  few 
matrix  multiplications  on  each  iteration.  For  consistency  in 
this  approach,  we  must  at  least  require: 


/{p*}j^{p*}j 


<<  1 


(163) 
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Apart  from  this  rule  of  thumb,  conditions  for  convergence  of 
(159)  remain  the  object  of  further  research.  At  present  we  can 
only  appeal  to  computational  experience  with  specific  cases  of 
which  a  detailed  report  will  be  given  separately.  As  (163)  would 
suggest,  (159)  is  usually  convergent  if  the  reciprocal  decorrela¬ 
tion  times  and  the  modal  frequency  differences  are  sufficiently 
large.  In  fact,  (159)  is  typically  linearly  convergent  for  all 
I  ^  0  if  the  modal  density  decreases  rapidly  enough  with  modal 
order  and  falls  below  some  bound  determined  by  the  other  con¬ 
ditions  of  the  problem. 

In  addition,  numerical  experience  or  inspection  of  (163)  give 
rise  to  the  following  conjecture.  With  the  partitioning; 


Pc 


Pp  e  C 


'CD 


2Nj^x2Njj 


(164) 


and  with  fixed,  there  is  an  large  enough  (and  comparable  to 
the  order  of  the  quasi-deterministic  system  which  is  not  directly 
affected  by  the  decorrelation  times)  that  the  sequence 


<Pktl>^m  =  [R  H  S  of  (162. a)]- 


m  =  2Np+l , .  .  .  ,2N^  ,  <  m 


^Pjc+l^Jl  =  [R  H  S  of  (162. b)]- 


(165) 


Z  =  2Njj+l, .  .  .  ,2N^ 


A 


is  convergent.  Here  denotes  (...)  evaluated  with 

fixed. 

If  true,  this  conjecture  suggests  the  following  scheme  for 
solution  of  (158)  when  is  large: 

(1)  To  start:  Set  =  0  in  (164)  and  determine 

by  the  methods  of  Section  3.6. 

(2)  With  Fp  as  determined  previously,  compute  and 

F^^^  by  use  of  (165). 

(3)  Given  F^,^  and  F^^^  as  computed  in  step  (2)  ,  calculate 
Pj^  from  the  upper  left  sub-block  of  (158)  by  either 
of  the  methods  of  Section  3.6.  Then  return  to 

step  (2)  unless  convergence  is  adequate. 

Consequently,  the  generally  convergent  but  more  elaborate 
Newton  method  or  perturbed  weighting  method  is  applied  to  a  rela¬ 
tively  low  order  stochastic  Riccati  equation.  At  the  same  time 
the  large  order  sub-blocks  F^p  and  F^^  are  handled  with  the  con¬ 
ditionally  convergent  but  much  faster  asymptotic  refinement 
method.  Thus,  although  general  conditions  for  the  convergence 
of  this  procedure  remain  a  subject  of  further  inquiry,  it  is  seen 
that  the  special  form  of  the  stochastic  Riccati  equation  admits 
the  possibility  of  efficient  numerical  treatment  for  systems  of 
rather  large  order. 
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4.  SUMMARY  AND  CONCLUSIONS 


This  work  has  addressed  two  of  the  principal  obstacles  facing 
the  application  of  modern  control  theory  to  structural  vibration 
suppression:  parameter  uncertainties  arising  from  the  intrinsic 
inaccuracies  of  structural  modelling  and  difficulties  in  the  form¬ 
ulation  of  optimal  control  laws  imposed  by  the  large  dimension  of 
the  system. 

At  the  outset,  the  need  for  a  design  method  which  takes  full 
advantage  of  the  peculiarities  of  the  system  to  be  controlled  in 
providing  a  statistical  treatment  of  a  priori  parameter  uncertain¬ 
ties  was  recognized.  In  addition,  it  was  emphasized  that  a  suit¬ 
able  method  must  employ  limited  data  on  parameter  statistics, 
thereby  eliminating  the  need  for  a  complete  probability  model  and 
reducing  the  number  of  required  elementary  measures  of  parameter 
variation  to  a  manageable  level. 

A  basic  inspiration  has  been  the  notion  that  by  consistent  use 
of  a  system  model  incorporating  limited  parameter  information  one 
can  so  arrange  matters  that  the  computation  required  for  formula¬ 
tion  of  an  optimal  control  policy  is  correspondingly  limited. 

Thus  it  was  hoped  that  the  obstacles  of  uncertainty  and  dimension¬ 
ality  could  both  be  circumvented  by  a  single  design  methodology. 

To  render  this  initial  development  tractable,  we  limited  con¬ 
sideration  to  full  state  feedback  regulation  in  the  presence  of 
uncertainties  in  the  open-loop  modal  frequencies.  In  addition, 
the  average  of  a  standard  quadratic  functional  over  both  disturb¬ 
ance  and  parameter  ensembles  was  chosen  as  the  performance  measure 
Within  these  restrictions,  considerable  progress  in  the  direction 
outlined  above  can  now  be  claimed. 

Chapter  2  set  forth  the  essential  ideas  of  a  new  approach. 

This  has  been  summarized  in  Section  2.4.  To  briefly  recapitulate, 
the  minimum  set  of  parameter  statistical  data  needed  to  preserve 


any  fidelity  in  the  overall  model  was  first  sought.  The  control¬ 
led  or  uncontrolled  system  was  seen  to  be  comprised  of  several  more 
or  less  distinct  qualitative  regimes  delimited  according  to  the 
values  of  the  "modal  decorrelation  times. "  These  parameters  con¬ 
stitute  new  measures  of  parameter  uncertainty  and  are  essential  to 
an  adequate  description  of  second  moment  response.  A  measure  of 
design  conservatism  is  then  achieved  by  constructing  a  full  proba¬ 
bility  model  of  frequency  uncertainties  which  is  unconstrained  save 
for  prescribed  values  of  the  decorrelation  times.  The  resulting 
white  parameter  model  allows  formulation  of  statistically  closed 
equations  determining  the  mean-square  optimal  control  law. 

At  the  close  of  Chapter  2,  the  problem  was  reduced  to  the  solu¬ 
tion  of  a  stochastic  Riccati  equation  which  assumes  the  form  of  the 
standard  Riccati  equation  when  all  parameters  are  deterministic. 

The  specific  properties  of  this  modified  Riccati  equation  were 
developed  in  Chapter  3.  First,  the  situation  with  regard  to  the 
existence  and  uniqueness  of  steady-state  control  laws  was  found 
to  be  entirely  analogous  to  that  of  deterministic  plant  regulation. 
Under  the  usual  stabilizability  and  detectability  restrictions, 
unique  steady  state  solutions  were  found  to  exist  for  all  levels 
of  uncertainty  in  the  modal  frequencies  -  indicating  the  absence 
of  an  uncertainty  threshold.  In  addition,  for  such  constant  gain 
controls,  strong  assurances  of  closed-loop  stochastic  stability 
were  found  in  Section  3.5. 

The  asymptotic  properties  of  steady  state  solutions  for  large 
uncertainties  and  high  order  modes  were  investigated  -  with 
significant  implications  for  the  dimensionality  issue.  In  the 
limit  as  all  uncertainties  increase  without  bound  the  solution  of 
the  stochastic  Riccati  equation  (in  the  eigen-basis  of  the  nominal 
system  matrix)  reduces  to  a  diagonal  matrix  whose  elements  are 
independent  of  modal  frequency  statistics.  This  asymptotic  solu¬ 
tion  gives  rise  to  a  velocity  feedback  control  law  which  is  stable 
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regardless  of  the  values  of  modal  frequencies  or  damping  ratios. 

If  all  uncertainty  levels  are  not  small  but  do  increase  with 
increasing  modal  order,  the  solution  approaches  the  asymptotic 
solution  for  high  order  modes  but  may  resemble  a  deterministic 
plant  control  for  very  low  order,  well-known  modes.  In  other 
words,  an  inherently  robust  velocity  feedback  control  automati¬ 
cally  emerges  for  the  poorly  known  high-order  modes. 

This  feature  greatly  reduces  dimensionality  problems  in  the 
solution  of  the  stochastic  Riccati  equation  as  shown  in  Section 
3.8.  The  convergent  numerical  procedures  set  forth  in  Section 
3.6  need  only  be  applied  to  the  "coherent"  portion  of  the  system 
which  consists  of  relatively  few  well-known  modes,  while  the 
velocity  feedback  control  for  the  high-order,  "incoherent"  portion 
can  be  determined  with  negligible  computational  burden.  Which 
modes  are  to  be  included  in  the  coherent  system  can  be  determined 
in  advance  of  any  burdensome  calculations  according  to  the  accu¬ 
racy  desired  in  the  determination  of  the  control  gain.  Finally, 
even  if  the  coherent  system  is  inconveniently  large,  the  computa¬ 
tional  schemes  advanced  in  Section  3.9  show  promise  as  efficient 
means  of  solution.  We  conclude  that  the  stochastic  Riccati  equa¬ 
tion,  which  results  from  inclusion  of  parameter  uncertainties  in 
the  fundamental  system  model,  is  amenable  to  satisfactory  numeri¬ 
cal  treatment  for  systems  of  arbitrary  order.  Detailed  numerical 
results  illustrating  the  above  features  will  be  given  separately. 

Obviously,  to  be  of  practical  value,  this  formulation  must  be 
extended  beyond  the  restrictive  assumptions  initially  adopted. 
First,  more  general  types  of  parameter  uncertainties  must  be 
accommodated.  Secondly  and  most  importantly,  the  assumption  of 
full  state  feedback  must  be  removed  and  the  theory  extended  to 
treat  fixed-order  dynamic  compensation  in  the  presence  of  parameter 
uncertainties.  These  generalizations  are  the  object  of  further 
investigation. 
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APPENDIX  1 


Proof  of  Theorem  1 

A.  First,  if  W(tj^,t2)  is  differentiable  almost  everywhere 
then  by  direct  substitution,  it  is  seen  that  (16)  satisfies: 


■^())(t,T)  =  [u  +  v(t)  -6k  (t)  ]  <|)  (t,  t) 
<l>  (T  ,T)  =  I 


(Al.l) 


almost  everywhere  and  is  the  transition  matrix  of  (11) .  If 
W(ti,t2)  is  almost  nowhere  differentiable  (see  References  32,  21), 
we  must  re-interpret  (Al.l)  as  * 


d<()(t,T)  =  [u  “6<  (t)  ]  <(i  ( t ,  t)  dt  +  dW(tH(t,T) 


(A1.2) 


where  dW(t)  is  the  differential  increment  of  W,  where  W  now  pos- 
esses  independent  increments.  (16)  may  be  shown  to  be  the  solu¬ 
tion  of  (A1.2)  in  the  sense  that 


/  d(j)  ( t ,  T ) 


=  [u  -6k  (t)  ]  (|>  (t ,  T  )  dt 
+  dW(t)((.(t,T) 


(Al. 3) 


holds  with  probability  1  for  all  a,b. 
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To  be  consistent  with  the  assumed  almost  everywhere  continu- 
ity  of  W(t),  we  must  take  W(t)  to  be  the  limit  of  a  band-limited 
process  as  the  total  power  approaches  infinity.  This  requires 
care  in  the  definition  of  the  stochastic  integral: 


t. 

G  =  /  dW(t)0(t,T) 

a 


(A1.4) 


appearing  in  (A1.3).  Following  Stratonovich ^ and  Wong  and 
Zakai^^^^,  the  proper  definition  of  (A1.4)  is  found  to  be 


G  4  Urn  2 


4>(t^,T)+(t>(t^^j^,T) 


(A1.5) 


as  MAX  (t^  , , -t  ,) -*0,  where  ft^,}  is  a  partition  of  [a,b].  With 

^  v  +  X  V  V 

this  interpretation,  (16)  is  seen  to  be  the  formal  solution  to 
(19).  Moreover,  <(i(t,T)  has  the  Ito  stochastic  differential  with 
Stratonovich  correction  term  of  the  form  given  in  Eq.  (20). 


Next,  the  series  of  (16a)  must  be  shown  to  be  convergent. 
Since  ►"(t)  is  bounded  and  W(T,t)  is  almost  everywhere  bounded, 
there  exist  finite  M^ ,  ”w  ^  0  such  that  for  all  t  t:(t^,tj^]. 


T  c [ t^ , 1 1  : 


6<(t)  I  I  <  M^,  I  Uo(t'T)  \  '  \ 


where  | |... | |  denotes  any  matrix  norm.  Then  it  may  be  shown: 
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A. _ t  , 


l\(t,T) 


P.l 


1 

<  kT 


(t-T) 


SO  that ; 


Utt.T)!  exp[M^M^(t-T)  ] 


Therefore,  the  series  (16a)  is  absolutely  convergent  and  ())(t,T) 
exists  and  is  bounded  almost  everywhere  for  all  finite  (t-t). 

B.  Since  <(t)  is  continuous  and  W(tj^,t2)  is  almost  everywhere 
continuous,  each  =  0,...“)  shares  this  property  as  is  evident 

from  (16b, c) .  Furthermore  (16a)  is  absolutely  convergent,  whence 
<l>(t,T)  is  itself  almost  everywhere  continuous  in  t  and  T.  This 
property  immediately  ensures  continuity  of  the  first  and  second 
moments  of  <ti(t,T).  To  show  differentiability  of  the  mean  of  <t)(t,T) 
we  first  average  (16c)  over  the  W  ensemble.  By  virtue  of  (16b)  it 
is  clear  that  each  element  of  E[4>j^(t,T)]  is  a  functional  of  Bk(Tjj^), 
Jl  =  l,...k  and  the  joint  characteristic  function  of  W(Tj^,t), 
W(T2,Ti) . • .W(T,Tk) .  But  this  characteristic  function  is  differen¬ 
tiable  in  all  its  arguments  since  increments  (15)  possess  first- 

/  OQ  \ 

order  absolute  moments. '  '  Moreover  the  arguments  corresponding 

to  W(Tj^,t),  W(t2  ,  Tj^)  .  .  .W (t  ,  Tj^)  are  proportional  to  (t-x^^)  , 

...(Tj^-t)  respectively.  Hence,  E[(J)j^(t,T)  ]  is  a  functional  of  con¬ 
tinuous  and  differentiable  functions  of  t,Tj^...Tj^,  and  consequently 
E(((»(t,T)]  is  continuous  and  differentiable.  An  analogous  argument 
shows  this  property  with  respect  to  second  moments  of  4)(t,T). 
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.J. 


C.  and  D.  Direct  substitution  of  (17)  and  (18)  shows  these 
expressions  to  be  a  formal  solution  of  (14) .  The  almost  every¬ 
where  boundedness  and  continuity  of  (j)(t,T)  shown  in  the  proof 
of  part  A  and  the  assumed  continuity  and  boundedness  of  6k 
ensure  the  same  properties  for  <^(t,T).  Consequently  expression 
(17)  exists,  is  bounded  and  continuous  almost  everywhere. 

Substituting  6  =  T-t  in  (18)  gives 

ij;(t,T)  =  (()“(T,t)[aj^  +  k”(t)R2k(t)  ]4>  (T,t)  (A1.6) 

Therefore,  since  >  0  and  R2  >  0,  >Jj(t,T)  is  positive  semi- 
definite  and  so  is  p(t). 

Interpreting  (14)  in  the  sense  of  (A1.2)  and  (A1.3),  and 
supposing  that  another  solution  distinct  from  (17)  exists,  we 
find  that  the  difference  obeys  the  homogeneous  form  of  (14b, c) . 
This  has  only  the  trivial  solution,  contrary  to  hypothesis  and 
it  follows  that  solution  (17)  is  unique. 

Finally,  from  (18),  each  element  of  E[ijj(t,T)]  ^  tJ;(t,T)  is 
given  in  terms  of  the  second  moments  of  il)(T,t).  From  part  B.  it 
follows  that  ^(t,T)  is  a  continuous  and  differentiable  function 
of  t .  □ 
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APPENDIX  2 


A 


Proof  of  Theorem  3 

Here  we  derive  the  result  of  Theorem  3.  First  consider  the 
discrete  problem.  The  entropy  may  be  defined  using  the  joint 
statistics  of  the  increments: 


^  ~  lf.../N  7  k  —  l/.../2n 


V  (A2.1) 


in  place  of  increments  (36).  With 


2n 


(X 


denoting  the  joint  density  of  increments  (A2.1),  where 
argument  corresponding  to  we  wish  to  determine  the 
which  maximizes 


the 

2n 

^N 


H  A  -/dX^...dX^"  pi  2n  iln  p  1  2n  (A2.2) 

subject  to  conditions  (a),  (b)  and  (c) ,  with  (35. a)  re-written  in 
terms  of  the  increments  (A2.1). 
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t4 _ 


The  constraints  imposed  by  the  available  information  are 

seen  to  be  independent  of  the  joint  statistics  of  any  two 

sequences  x^-*  )!,  =  1,...N  and  xl’  ^  1,...N  for  k  /  j.  More- 

^  k  ^  i 

over,  H  is  maximized  for  Xo  and  xi  >  t  ^  independent.  We 

1  2 

may  conclude  at  once  that 


2n 

pi  2n  =  n 
Xi-"Xn  k=l 


H 


=  E  H" 
k 


t  -/dX^. 


•dX^ 


k  «,n  p^k 


Xi* • -X^ 


y  (A2.3) 


k 

and  it  remains  to  determine  each  p  k  k  to  maximize  H 
,  Xi*--Xn 

separately. 

Consider  k  fixed,  and  drop  the  k  superscript.  Ws  may  write: 


P 


X 


1* 


•  x 


N 


r 


L  (Xt 
Xfj  Xi  2 


X„|Xi)/p^, 


( X2 .  .  .  Xjj) 
Xjj  ^ 


L  (A2.4) 


introducing  r  as  a  new  independent  function  embracing  the  statis¬ 
tical  dependence  between  Xj^  and  (X2  » •  •  • 'Xjj)  •  Clearly  F  must 
satisfy  the  constraints 
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/dx,  r  p  =1 

/dX- . . . dx„  r  p  =  1 

2  N  ^X2---Xn 


(A2,5) 


Now  we  may  append  (A2.5)  to  the  original  problem  and  seek  F, 
p  and  p  to  maximize  H  .  If,  for  the  moment,  condition  (a) 

Xi 

is  removed  from  consideration,  it  is  easily  seen  that  F  =  1  maxi¬ 
mizes  the  entropy.  With  this  value,  (35)  may  be  written: 


1 

N 


[1  + 


N 

2 

m=2 


‘2m' 


■]  = 


Tk/T 


(A2. 6) 


where 


A,  ^  /dX,  p  e^^'k^l 
1  “  1  Xi 

_  m 

4  /dXj.-.dX^  e  ’‘.•=2 


(A2.7) 


and  Tj^/T  <  1  by  (35b). 

Since  Aj^  and  A2j^^  are  characteristic  functions 


1Ai|2,  |A2„|"  £  (0,  1) 
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so  that  there  always  exists  a  choice  of  p  and  p  such 

Xi 

that  (A2.6)  is  satisfied.  Consequently,  is  maximized  under 
all  constraints  by  F  =  1,  and: 


But,  by  the  stationary  condition  (b) : 


Py  y  (^2'  '  ’  Py  y  (X  X  i 


Xl  ^  X2-**Xn_i  ^ 


Repetition  of  this  argument  (N-3)  times  shows  that 


P 


X 


1* 


•  •  X 


N 


N 

n 

m=l 


Since  this  choice  maximizes  the  entropy,  we  conclude  that 
(A2.2)  is  maximized  under  conditions  (35),  (b)  and  (c)  if  the 
increments  Xo  are  all  mutually  statistically  independent  and  if 

Jr 

for  each  k,  I  =  1,...N  are  identically  distributed.  This 

result  holds  for  all  N,  T  >  Tj^,  k  =  1...2n. 
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Beyond  this  point,  further  conditions  for  entropy  maximiza¬ 
tion  need  not  be  considered,  because  in  the  limits  and  Tt<», 

(35)  uniquely  determines  the  statistical  structure  of  the  6j^{t). 

Considering  the  limit  N+“,  it  is  found  that  the 
t  e(0,  T)  are  mutually  independent  and  possess  independent 
increments.  Any  increment  may  be  expressed  as  the 

limit  of  the  sum  of  sub-increments  defined  over  equally  spaced, 
disjoint  intervals.  The  sub- increments  are  independent  and  by 
conditions  (b)  and  (c)  are  stationary  and  possess  finite,  non¬ 
zero  variances.  Thus,  the  conditions  of  the  Lyapunov  central 
limit  theorem hold  and: 


6j,(t^,t2)  =  N(0,Oj^(tj^,t2)) 

k  ”  1,..., 2n 


(A2.8) 


where  denotes  the  variance.  Furthermore,  that  '5j^(t)  is  zero 
mean  and  its  increments  are  stationary  and  independent  suffice 
to  imply 


^)c^^l^^2^  —  I  ^2  ”^1 


(A2.9) 


where  the  are  positive  constants. 

With  (A2.8)  and  (A2.9),  (35)  yields: 


1  -u)„  .  B,  T  1  , 

_ I  _  i  N(k)  k  _  A  1 

32  X  e  -  =  32 - 

‘^N(k)®k  ‘^N(k)^k 
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in  the  limit  N  +  «>.  This  determines  uniquely.  In  the  limit 
Tt» : 


u). 


/  k  —  I/,,,/ 2n 


N(k) 


(A2, 10) 


and  (39)  follows. 
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APPENDIX  3 


A 


Proof  of  Corollary  2 

First,  it  is  advantageous  to  rearrange  the  state  vector  so 
that  modal  coordinates  are  the  first  n  states  and  modal  velocities 
are  the  second  n  states.  In  other  words,  perform  the  coordinate 
transformations 


E 

E 

o 

o 

X  = 

E 

X'  ,  ^  = 

E 

e 

e 

(A3.1) 


th  T  th 

where  the  m  row  (m  =  l,...,n)  of  E  is  ei;  ,  while  the  m  row 

o  2m— 1 

of  E^  is  el„  where  e,  denotes  the  standard  unit  basis  vector, 
e  2m  1 

Under  this  transformation,  the  various  matrices  appearing  in 
(48)  become: 


V  o 

O  V* 


V 


A 


diag  ,  iw2“n2‘^2 


•)  ^ 


I 


J  o 

o  J 

L 


/  U  ^  diag(u)2^l2,^2^^^***) 


a 


1 


’12 


'12 


o 


2 


C  -C 
-C  C 


>  (A3. 2) 
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where  C  is  real.  Consistent  with  the  above  partitionings,  the 
matrix  of  "p  under  transformation  (A3.1)  may  be  written 


P 


*^11  ^12 
^12  ^22 


(A3. 3) 


so  that: 


J{p} 


0  '^^^22^ 


(A3. 4) 


Using  (A3. 2),  (A3. 3)  and  (A3. 4)  we  may  expand  (48)  into  its  four 
sub-blocks.  Manipulation  of  these  relations  yields: 


•  H 

■z-^  =  (v-*5J)Zj^  +  Zj^(v-J5J)  ~Zj^Cp2^j^-P22^  ^2  ^1*"  ^2 


^1^  ^12  ^12^  ^1  "^l^^  ^2^Pll  Pll^  ^2 


^2^  Pl2"Pll‘^  ^1 


y  (A3. 5) 


{ 
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r 


-22  =  (V*  -h3)Z2  +  ^2^^ 


-Z2C  Pii-PiiC  Z2  +  Z2C  Z2  +  Z2C  PJ2 

+zjc  Pj_j_  +  pj_2C  Z2  -z“c  Z2  -zjc  PJ2 


+  pj^j^C  z-^  -Z^C  Zj. 


P12C  z”  +  z^C  z^ 


MA3.6) 


Z2(ti)  =  0 


where : 


z  A  n * 

=  ^12  ^12 


Zo  ^  P 


-  n* 


11  "^22 


(A3. 7) 


From  the  homogeneity  of  (A3. 5)  and  (A3. 6)  in  Zj^  and  Z2  we  conclude: 


Zj^  —  Z2  —  0  ;  te[ t^ f 1 2^ } 


(A3. 8) 


so  that  any  hermitian  solution  of  (48)  assumes  the  form: 


Pll 

P12 

p  = 

_Pl2 

Pll 

Pll  = 

Pll  / 

II 

CM 

rH 

Q. 

y  (A3. 9) 


b .  /  c . 
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Under  transformation  (A3.1),  the  definition  (8)  produces; 


$  = 


I  I 
ifi  -in 


(A3. 10) 


n  ^  diag  (a)j^  ,012  / .  .  . 


so  that 


fT  A  .-1 

P  £  4i  p4> 


=  H 


Re[Pij^+Pl2] 


-1 


-n  ^Im[p2^j^+p2^2l  ^  ^Re[Pj^]^-P ^^2^ ^  ^ 


{A3. 11) 


This  shows  that  P  is  real  and,  by  virtue  of  (A3.9b,c)  is 
symmetric.  □ 
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APPENDIX  4 


Proof  of  Lenuna  1 
— 

A.  With  )  detectable,  the  rigid  body  modes  are  neces¬ 

sarily  reconstructible.  Furthermore  since  dcimping  in  all  elastic 
modes  is  assumed  non-zero,  only  states  corresponding  to  elastic 
modes  are  contained  in  the  unreconstructible  subspace.  Now 

is  diagonal  and  its  non-zero  elements  distinct  since  the 
w,  ,  k  =  l,...,n  were  assumed  distinct  in  Section  3.1.  Thus  the 
unreconstructible  subspace  is  spanned  by  a  set  of  unit  basis 
vectors  (in  the  eigen-basis  of  A)  corresponding  to  the  distinct 
unreconstructible  poles  of  With  an  analogous  result  for 

the  reconstructible  subspace,  it  follows  that  the  system  equa¬ 
tions  may  be  put  into  a  reconstructibility  canonical  form  merely 
by  rearrangement  of  the  state  vector  as  in  (79) .  The  forms  given 

by  (80)  for  and  o,  then  follow  at  once, 
ml 


B.  First  we  show  that  the  solution  to  (64)  exists.  The  right 
side  of  (64a)  is  analytic  in  ^j^j,V(k,j),  whence  p(t)  is  continu¬ 
ously  differentiable  at  t  when  p(t)  exists.  Hence  p(t)  also 

(48) 

exists  on  [t  -e,T]  for  some  e  >  0  by  analytic  continuation 
Then  since  ^(t^^)  is  given,  p(t)  exists  on  (-“>,tj^]. 

Next,  let  q  e  C  be  defined  on  (-“rt^^]  by 


•A 


so  that: 


q“(t)p  {t)q(t) 


tl  H 

;  dT  q“(T)  [I{p(T)  } 
t 


+  Oj^lqCT) 


and  it  thus  follows; 


q“(t)p(t)q(t)  >  0  ,  V  t  ^  tj^ 


Note  that  q(t)  can  be  selected  as  any  vector  in  because 

q(t)  =  'l>(t  and  state  transition  matrix  it  is  nonsingular. 

Thus  p(t)  is  positive  semi-definite  on  (-«>,t^] - 

Finally  we  show  that  F  reduces  to  the  form  given  by  (81)  and 
(82).  Partition  ^  in  a  manner  consistent  with  (79): 


Then  with  (80)  and  (83a)  the  (1,2)  and  (2,2)  sub-blocks  of  (64) 
may  be  written: 
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■^12  “  *^12  ‘^12  ^mu  "^f°2  ^12  ®  ^2  ®u  ^2^ 

-»12l«u  e"  ?12  +  B„  8“  JjJ 


■^2  ‘^mu^2  ■■  P2^mu  •^u^‘^2^  ~^12^^2  ^12  ®  ^2  ®u  ^2^  ^(A4.2) 


-P2(B„  R-^  b"  p,2  +  6„  r;^  bJ  Pj] 


''l2^h^  =  ^2^4^  =  ° 


It  is  seen  from  the  above  that  'p^2  ^2  continuously 
differentiable  at  t  =  t^^  and,  moreover,  all  derivatives  vanish. 
Thus  by  analytic  continuation: 


Pl2 


=  P2 


0  ;  V  t  <  tj^ 


and  (81)  is  justified.  Substitution  of  (81)  into  (64)  produces 
(82)  as  the  only  non-zero  sub-block.  □ 
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APPENDIX  5 


Proof  of  Theorem  8 

First,  note  that  by  virtue  of  Leirana  2  the  steady  state  solu¬ 
tion,  "p,  of  (64)  is  given  by 


0 

0 


Next 

definite 


define  Q(t)  by  (69)  with  v 
matrix.  Partitioning  Q(t) 


=  0  and  any  positive  semi- 
consistently  with  (79) ; 


(A5.1) 


and  substituting  this  into  (69)  along  with  (80)  and  (83)  ,  one 
obtains : 


^^m  ■‘^2^^®12  ^12  ^mu 

*  02  “Su  *  ^u<02' 

-B„  R-^  ?  A 


"  +  J{Q} 

-Q  A  6  Rj^  S" 

-qJj  a  5  R]^ 


(A5. 2) 
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--  H  ^ 

By  assumption,  -a2A)Q  +  Q(Uj^  -^2^)  +  -^(Q)I  is 

asymptotically  stable.  Thus,  from  (A5.2a),  Q  converges  expo¬ 
nentially  to  zero  and  there  exist  real  a  ,6  >0  such  that: 

o  o 


o  -a^(t-t  ) 
,  <  6^  e  o'  o' 
2.  O 


{A5. 3) 


Note  that  the  contribution  of  T{q]  to  is  of  the  form  of 
a  diagonal  matrix,  each  diagonal  element  being  either  zero  or 
equal  to  one  of  the  diagonal  elements  of  r.  Consequently: 


-  '2{u„  -02A)“  +  /(Q}]) 


where  denotes  the  real  part  of  the  eigenvalue.  Since, 

by  assumption  the  right  side  above  is  negative, 

A  A  A  A  A  ^  „ 

A^[(ir  -o_A)Q  +  Q(iJ  -o_A)"]  is  asymptotically  stable;  and  this 
implies  that  -a^A  is  asymptotically  stable. 

By  the  detectability  assumption  of  Lemma  1  and  the  condition 

_ U 

of  non-zero  damping  on  all  elastic  modes,  is  also  asymptoti¬ 

cally  stable,  whence  (A5.2b)  has  the  unique  solution: 


(Um-02A)(t-to) 
Qi2  =  e  Q 


12^^o' 


-  f^d-x 


(lJ„-o,A)  (t-T)  - 
e  ^  ^  Q  A 


6  R 


-1 


-H  , 


u 


124 


-  Am 


The  stability  of  ^mu  exponential  convergence 

to  zero  of  the  first  term  above.  Moreover,  the  bound  (A5.3)  and 

A.  /S 

the  stability  of  (y^  ^mu  exponential  convergence 

for  the  second  term  as  well.  Therefore,  there  exist  real  and 
positive  01^,  such  that: 

(t-t  ) 

1 \Q^2^  I2  ^  ^1  ®  (A5.4) 


_  _tJ 

Next  we  show  that  A,.,-[y„,0.,  +  Q-y  +  I  {Q-,}]  is  asymptoti- 

Q2  mu  2  2  mu  u  2  ' 

cally  stable  for  all  I  >0.  Consider 
^  u 


(AS. 5) 


for  any  hernitian  Q„  >0.  Recalling  that  y  =  u  -*5/  ,  we  have 
^  2o  mu  u  u 

from  ( A5. 5)  : 


■=2  -  ‘•r.u02  ^ 


Q2<^o>  =  Q 


20 


^  tr  Q2  =  2  tr [ (Re  u^)Q2l  <  0 


where  the  equality  on  the  right  holds  only  if  1 IQ2I I2  “  More¬ 

over,  as  (tr  ©2 J  '  0  for  t  clt^,®),  tr  Q2  is  seen  to  converge 
exponentially  to  zero  as  t  +  “>.  This  implies  exponential  conver¬ 
gence  for  fQ2^  also  since  |l{Q2’ll2  ^2'  (A5.5),  Q2(t) 

satisfies 
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.. 


u  (t-t  )  (t-t  ) 

^  ^mu  o'  ^mu  o' 

Qlit)  =  e  02^6 


+  /^  dT  e 
t 


(t-x) 


Iu{Q2>  e 


-H  X 

(t-x) 


Stability  of  and  exponential  convergence  of  {Q2)  imply 

exponential  convergence  of  Q2(t)  as  defined  by  (A5.5).  There¬ 
fore  A._[u^  0-  +  is  asymptotically  stable. 

Finally,  from  (^5.2c)  this  property  together  with  the  norm 
bound  (A5.4)  and  implies  asymptotic  convergence  of  IIQ2II2' 

Combining  the  above  results,  one  sees  from  (A5.1)  that 


LIM  Q  =  0 
t-t 


so  that  A„I(u  -a,p)Q  +  Q ( -op)”  +  J{Q}]  is  asymptotically 
Q  m  2  m2 

stable.  n 
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APPENDIX  6 


Proof  of  Lemma  6 

To  prove  the  first  statement,  consider; 


0  =  u^Y  +  Yu  +  -Y  02  Y 


-  0,  (u,  s^)  completely  reconstructible 


(A6. 1) 


The  conditions  on  s^^  and  the  assumed  stabilizability  of  (u,3) 
ensure,  by  Lemma  5  that  (A6.1)  has  a  unique  positive  definite 
solution  and  that  (u  “O2Y)  is  asymptotically  stable. 

“  1  ^  H 

Now,  with  the  control  gain  R-  B  Y,  the  reconstructible 

A 

State, C,  is  given  by: 


^  =  (U  -O2Y)  +  Z  (A6.2) 


where  is  zero  except  for  the  (k,k)  element  which  equals 
ilmlT.  .  Clearly  B,  is  skew-hermitian.  This  fact,  together  with 

K  Kaa  /  a  \ 

the  stability  of  p  -u^Y  imply'  that  the  null  solution  of 
(A6.2)  is  P  mean  exponentially  stable  for  any  P  (in  particular, 
P  =  2) .  Since  is  the  hermitian  transpose  of  the  equivalent 
coefficient  matrix  of  the  second  moment  matrix  of  (A6.2) ,  ' ,  is 
asymptotically  stable. 
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In  view  of  the  stability  properties  established  above  for 
(A6.2)  we  can  show  the  result  for  the  case  in  which  the  system 
lacks  rigid  body  modes  (i.e.  ,  vT  =  y  )  ,  by  demonstrating  the 
asymptotic  stability  of  u  ^  positive  and  diagonal. 

Thus,  consider: 


I  =  (y  -a^Y)  ? 

(A6. 3) 

=  ^0 

for  any  ^  0  positive,  diagonal  Y.  Defining 

^  k  C  (A6.4) 


we  obtain  by  (A6.3) : 

(Y*^C)“[Rey  -Y^a^Y^]  (Y^^)  (A6.5) 

dt  2 


The  right  side  is  negative  for  all  C  0  since  damping  is 
assumed  non- zero  for  all  elastic  modes.  Therefore  is  a 
Lyapunov  function  for  (A6.3)  and  y  ~02Y  is  asymptotically  stable. 
By  the  argument  following  (A6.2),  the  asymptotic  stability  of 
is  established.  □ 
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APPENDIX  7 


Proof  of  Theorem  11 

First,  (-JjIP  +  iCp}  -J3PJ)  is  shown  to  be  negative  definite. 
Let  W  be  a  real  diagonal  matrix  of  zero-mean,  statistically 
independent,  Gaussian  random  variables  such  that: 


E  [w: 


kk^ 


=  J 


kk 


{A7. 1) 


Then,  defining: 


r 


iWe—  -iWe 
p  e 


e  >  0 


(A7.2) 


one  obtains  the  identity: 


-^IP  +  l{p}  -Ho I  =  -LIM  A  E[ri  (A7.  3) 

e  +  O 

Employing  the  reconstructibility  decomposition  afforded  by 
Lemma  1,  T  assumes  the  form 


A 


where  A  is  the  positive  definite  solution  to  (88).  Letting 
W  ^  A^WA”**  and  C  =  A^C  for  any  C  ^  0,  we  have: 


(A7.5) 


Since : 


S  I  5 


where  the  last  line  follows  because  iW  is  skew-hermitian ,  (A7.5) 
yields: 


>  0  ,  V  C  7^  0 


Thus  r  and  therefore  T  is  positive  semi-definite  for  each  reali¬ 
zation  of  W  over  its  statistical  ensemble.  Clearly  E[r]  is  also 

2 

positive  semi-definite  and  is  of  order  e  by  definition.  Conse¬ 
quently,  the  right  side  of  (A7.3)  exists  as  a  negative  semi- 
definite  matrix.  This  suffices  to  show: 

S  <  (A7.6) 


i 


Next,  rearranging  (102) ,  p  is  seen  to  be  a  positive  semi- 
definite  solution  of 


0=u^P+PU-pa2p+S  {A7.7) 

Because  -  0  and  ^  s  0,  (IT^  ^  +  p  y)  is  negative  semi- 

definite.  Consequently: 


S 


-(u”  p  +  p  u) 


P  02  P  >  0 


{A7. 8) 


Finally,  we  show  that  detectability  is  preserved  despite  the 
presence  of  the  stochastic  terms  in  (107) .  Denoting 
{-Hi'p  +  l{p}  -Hpl)  by  a  and  partitioning  in  accordance  with  (59), 
one  obtains 


(A7.9) 


a  = 


0 

H 

a 

L  re 


by  virtue  of  (63a)  .  Denoting  the  "r"  sub-bloc)c  of  a,  by  o^  ,  it 

is  ^ 

follows  that  (Oj^  +  a)  assvimes  the  form 


(Oj^  + 


a)^  =  (a 


Ir 


s) 


(A7. 10) 


i .  e 


the  first  2n^  columns  of 


{o^  + 


a)  ^  and  of 


o 


*5 

1 


•  t 


are  identical. 


A 


Now,  with  of  the  form  (60) ,  the  necessary  and  sufficient 
condition  for  complete  reconstructibility  of  the  rigid  body  modes 
is  (see  p.45  of  Reference  (44)). 


rank[vr^,(a^  +  a)^^]  =  n^. 


(A7.11) 


where  (o,  +  a) ^  denotes  the  "r"  sub-block  of  (o  +  a)^,  and  is 
merely  by  (A7.10).  Thus  (A7.11),  which  becomes: 


— H 

rank , 


“Si  ■  "r 


(A7.12) 


—  h 

is  satisfied  since  (u,  o?)  was  assumed  detectable.  Therefore 

_  k 

(M,  S^)  IS  detectable.  □ 


APPENDIX  8 


Proof  of  Theorem  17 

The  boundedness  conditions  on  B  and  R2  may  be  stated  as: 


’’kj  ‘  ^ 

a. 

^2kk  ^  ^ 

b. 

n  <  rij^  <  ^ 

c. 

>  (A8.1) 


where  b,  r,  n  and  n  are  positive  and  the  b,  .  are  non-zero  ele- 
ments  of  B  as  given  by  (3).  From  (3),  (8)  and  (A8.1a,b)  it 
follows: 


2kk 


E 

4r  ^ 


(A8. 2) 


N(k) 


This,  together  with  condition  (2)  of  the  Theorem  and  (123)  yield: 


l<P\jl 


<  ‘*^l'^N(k)  ‘^N(j) 

^k‘^N(k)  ^j‘^N(j) 


(A8. 3) 
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with  use  of  (A8.2)  and  (A8.3),  manipulation  of  (124a)  produces: 


°i*k 


Similarly,  (124b)  becomes: 


.4-2-2  r-2^2 

"  r  °1  ‘^N(lc)  ^ 


(A8.4) 


(A8. 5) 


where 


f  ^  2  _ _ _ t _ - _ 

^  ^1<“n(1c)  ^ji“n(ji) 


(A8. 6) 


In  view  of  conditions  (1)  and  (3)  ,  possesses  an  upper 

bound  proportional  to  the  Riemann  zeta  function,  C(s) ,  with 

s  >1.  Therefore  f.  is  convergent  as  nt»  and,  moreover,  both 
LIM  ___  _ _ _ _ ^ _ .. _ 4 _ 


fk  and 


nt®  monotone  decreasing  with  increasing  k.  Con¬ 


sequently,  there  exists  an  N  sufficiently  large  that: 


{Ol}k  ^  0 


r-2 
On  b 

1  -  -  f  >  0 

■Q^‘^N(k) 


(A8.7) 
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for  all  k  >  2N,  and  we  may  restrict  attention,  in  the  following 
to  the  behavior  of  and  for  N  >  N. 

X  C 

Condition  (A8.7a)  permits  the  use  of  part  (b)  of  Theorem  (15) 
in  the  evaluation  of  an  upper  bound  to  I  But  first, 

bounds  must  be  determined  for  Note  that  (122)  can  be 

written  in  the  form; 


’Ikk  ,1 


{p*},.  = 

^  I  Reuj^  I  * 


A  ''lkk°2kk  ^  Q 
=  (Re  V2  ■ 


b. 


>  (A8.8) 


Obviously ; 


{p*}k  5  2 


1  ^Ikk 


I  Reuj^  1 


,  V  k,  X  >  0 


By  virtue  of  (A8.2),  (A8.1c)  and  condition  (2)  of  the  Theorem, 

X  decreases  with  k  increasing.  Consequently,  as  (A8.8a)  reveals, 
for  given  N  there  is  an  M  >  1  such  that; 


{p*} 


°lkk 


V  k  >  2N 


In  the  following  discussion,  N  is  considered  fixed  with  N  <  N  so 
that  M  is  fixed.  With  this  proviso,  we  have; 
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r 


-1  - 

2M  -  ‘^N(k) 


1  ^l^^Ndc) 

2  n 


V  k  >  2N 


(A8. 9) 


Using  (A8.1),  (A8.4),  (A8.5),  (A8.9)  and  conditions  (1) 
through  (3)  of  the  Theorem,  (126)  may  be  made  to  assume  the 
form: 


{p}k 

— — 


<  n  M  -  r-2 

^  - - a,  b 

n  r  1 


^  ‘^N(k) 


1  - 


r-2 

Oj^b 


‘^N(k) 


(A8.10) 


Since  fj^  is  monotone  decreasing  with  k  and  in  view  of 
{A8.7b),  the  right  side  of(A8.10)  is  a  positive  monotone  de¬ 
creasing  function  of  k  for  k  >  2N.  Thus,  given  e  >  0,  there  is 

an  N  such  that; 
c 


s  e 


k  >  N 

c 


(A8.11) 


In  particular,  let  correspond  to  £  <  1  in  {A8.11). 
Then : 


{p}jj  ^  {p*}j^(l  -g)  >  0  ;  k  2 


(A8. 12) 
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Using  this,  (A8.3)  yields: 


]<^>  .1  ,  /  2Mn  ‘^N(k)‘^(j) 

V{p}j(l-S)£^/  VN(k)  ^j“N(j) 


k  S:  N 


cl 


(A8. 13) 


The  right  side  of  the  above  inequality  is  a  positive,  mono¬ 
tone  decreasing  function  of  k  for  k  ^  Therefore,  from 

(A8.11)  and  (A8.13),  it  is  possible  to  determine  a  ^ 
such  that  both: 


Hp>„  s  c  ,  V  k  ^ 

and 

|<p>,^j|  s  e  ;  V  j  ,  V  k  > 


(A8. 14) 


These  relations  suffice  to  show  (152)  directly.  □ 
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